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ABSTRACT 

Context. The spatial distribution and variability of Fe-Ko- emission from molecular clouds in the Galactic Centre region may provide 
an important key to the understanding of the recent history of Sgr A*. A very plausible interpretation is that this variability represents 
an echo in the reflected radiation from the clouds of a past episode of high activity in Sgr A*. 

Aims. We examine the temporal and spectral properties of nine Fe-Ka bright molecular clouds within about 30 pc of Sgr A*, in order 
to understand and constrain the primary energising source of the Fe fluorescence. 

Mefhods. The variability of the Fe-Ka- line at 6.4-keV was investigated by spectrally fitting the data derived from the EPIC MOS 
cameras, after subtracting a modelled background. We have also studied the reflection imprints in time-averaged pn-spectra of each 
cloud, by measuring the equivalent width (EW) of the 6.4 keV line and the optical depth of the Fe-K absorption edge at 7. 1 keV. 
Results. Significant Fe-Ko' variability was detected, with a spatial and temporal pattern consistent with that reported in previous 
studies. The main breakthrough that sets our paper apart from earlier contributions on this topic is the direct measurement of the 
column density and the Fe abundance of the MCs in our sample. We used the EW measurements to infer the average Fe abundance 
within the clouds to be 1.6±0. 1 times solar. The cloud column densities derived from the spectral analysis were typically of the order 
of 10"^ cm"-, which is significantly higher than previous estimates. This in turn has a significant impact on the inferred geometry and 
time delays within the cloud system. 

Conclusions. Past X-ray activity of Sgr A* is the most likely source of ionisation within the molecular clouds in the innermost 30 
pc of the Galaxy. In this scenario, the X-ray luminosity required to excite these reflection nebulae is of the order of lO^'-lO^* erg s"', 
significantly lower than that estimated for the Sgr B2 molecular cloud. Moreover, the inferred Sgr A* lightcurve over the past 150 
years shows a long-term downwards trend punctuated by occasional counter-trend brightening episodes of at least 5 years duration. 
Finally, we found that contributions to the Fe fluorescence by X-ray transient binaries and cosmic-ray bombardment are very likely, 
and suggest possible ways to study this latter phenomenon in the near future. 
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' 1. Introduction 

I The Galactic Centre (hereafter GC) region is a unique environ- 
ment within the local Universe, which provides many tests of 
[ our understanding of fundamental issues in astrophysics. The re- 

■ gion hosts the nearest Super Ma ssive Black Hole (SMBH), Sgr 

■ A*, w ith a mass of 4xlO^M0 t Schodel et al.ll2002t iGhez et alJ 
; I2008h . It is also a region in which high energy phenomena 

■ abound. For example, in X-rays the Sgr B2 and Sgr C molecu- 
lar complexes, located at projected distances of 90 pc and 70 pc 
from Sgr A*, shine brightly through 6.4-keV Fe-Ka Une emis- 
sion ("Kovam a et al.ll99 6VMurakami et al.'200tf,'Nakaiima et al.' 
[2009), consistent with the prediction of Sunvaev et al. ( 1993). 

The physical mechanism responsible for the Fe-Ka emis- 
sion from these molecular clouds near the GC is the fluorescence 
of cold, neutral or near-neutral matter irradiated by high energy 
particles or X-ray photons. So far several hypotheses have been 
proposed as to the nature and origin of the primary source of 
this irradiation. The p ossibilities include: the X-r ay reflection 
nebulae model (XRN, ISunvaev & Churazovlll998l) . heating by 
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low-energy c osmic-rays (CR) (lYusef-Zadeh et al ] l2002h . shock 
mechanisms (Yusef-Zadeh et al. 1997), and electron bombard- 
ments (P redehl et al. 2003) . A recent suggestion is that subrel- 
ativistic protons can be created via accretion of stellar debris 
onto the central black hole, thus explaining both the observed 
X-ray continu um and the 6.4 keV line emission from the GC 
(iDogiel et al.l 12009). Other particle-like candi dates ca n be CR 
electrons originating in supernova events (. Valinia et al.ll2b00) . 

If the fluorescence observed in molecular complexes is the 
result of irradiation by X-ray photons, a luminous localised 
source of X-rays must be invoked to power the observed emis- 
sion, since the diffuse hot plasma which permeates the central 
regions of the Galactic plane produces a factor ten less photons 
than is required to account for the observed Fe-Ka flux. It has 
been estimated that, depending on its position relatively to the 
Sgr B2 and Sgr C molecular clouds, this source should have a 
2-10 keV luminosity of the order of 10^'' erg s '. Presently no 
persistent sources in the GC region have such a high luminos- 
ity; however a past transient outburst in a source, which has now 
returned to a relatively quiescent state, might well match the re- 
quirement. Sgr A* itself is arguably the best candidate, since 
although its activity is currently rather weak (the brightest flare 
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measured with X MM-Newton has a luminosity of the order of a 
few 10^^ erg s Porauet et all2003h . the SMBH could have un- 
dergone a period of high-state activity in the past ( Kovama et al] 
Il996 ). In fact, the projected light travel time to Sgr B2 implies 
an outb urst roughly 75-150 years ago, according to a recent es- 
timate (lTerrieretal.ll2010h ." Since there are a number of fila- 
ments nearer to Sgr A* which emit at 6.4-keV, the same argu- 
ment would require Sgr A* to also have flared more recently, 
say within the last 100 years, depending on the relative position 
of the clouds along the line of sight. 

One method of distinguishing between the point source and 
the particle hypothesis is to investigate the lightcurve of the fila- 
ments which emit the 6.4-keV line. If we are dealing with a tran- 
sient source as an engine of the fluorescent emission, one might 
expect the temporal evolution of the observed line flux to follow 
the light curve of the outburst, subject to blurring arising from 
the spread in light-travel delays inherent in the source to cloud 
geometry. Such effects have been found in the Sgr B2 molecular 
cloud where the bri ghtest peak of the Fe-Ka emission observed 
in 1995 by ASCA kovama et al][T996t [Murakami et al.ll2000h 
maintained the same brightness level when remeasured in 2000 
with Chandra dMurakami et al.l l200D but then, some five years 
later when observed by Suzaku, showed a marked decline to ap- 
proxi mately half t he peak value (Kovama et al.ll2008t llnui et all 
12009 ). Similarly, Nakaiima et al.l (i2009l) discovered variability 
of the Fe-Ka flux from Sg r C on the bas is of recent Sukaku ob- 
servations. Moreover. iMuno et al.l (l2007h have reported the evo- 
lution in intensity and morphology of the 4-8 keV continuum 
emission in two filamentary regions located close to Sgr A*. The 
latter changes occurred on parsec scales (in projection), which 
requires a brightening/fading of the illuminating source over a 
2-3 year period, with a n inferred 2-8 keV luminosity of at least 
10" erg s Recentlv. IPonti et all ( 1201 Ol) showed that the 6.4 
keV line flux from the molecular filaments with 15 arcmin of 
Sgr A* exhibit a complex pattern of variability. If these molec- 
ular clouds/knots have a particular distribution along the line of 
sight, then it is possible to argue that they were all energised by 
the same ou tburst on Sgr A*, consistent with the XRN scenario 
dPontietaLl UoiO). However, given the observational complexi- 
ties, it is very poss ible that this is not the complete story; indeed, 
very recentlv, Ca pelli et al.l ( [201 lb) studied the Fe-Kc line emis- 
sion from the MCs in the Arches cluster region (about 20 pc in 
projection from Sgr A*), showing that the XRN/Sgr A* scenario 
can hardly describe the spectral and temporal properties of those 
clouds. 

The above results open once more the question whether Sgr 
A* has exhibited AGN activity in the past and, if so, what is the 
exact nature of that activity. It is in this context that we reassess 
in this paper, the morphology, variability and spectral properties 
of the 6.4-keV emitting clouds within 15'of Sgr A*, using the ex- 
tensive set of XMM-Newton observations targeted at this region. 
Our goal is to investigate both the past role of Sgr A* (or other 
transient sources) in illuminating the GC molecular clouds and 
also to seek evidence for a contribution from alternative mecha- 
nisms, such as CR bombardment, in the excitation of the Fe flu- 
orescence which characterises the GC environment. Throughout 
this work, the distanc e to the GC has been taken to be 8 kpc 
dGiUessen et al.ll2009h . 

2. Observations and Data Reduction 

We checked in the HEASARC archive for XMM-Newton 
observations of the GC region and selected observations 
with a nominal pointing position within 15 'from Sgr A* 
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Fig. 1. Soft proton flai'e filtering of the MOSl dataset from OB SID 
0202670701. Upper panel: 2.5-12 keV count rate histogram. The blue 
lines mark the region used in the Gaussian fit, the green line represents 
the best fit Gaussian and the red lines show the bounds used to filter the 
data. Mid panel: 2.5-8.5 keV lightcurve of the IN FOV region. Lower 
panel: 2.5-8.5 keV lightcurve of the corner data. In both the Mid and the 
Lower panels, the points coloured green in the light-curves correspond 
to the selected GTI intervals (see text). 



(R. A.= 17''45'"40.045 ', DEC.=-29°0'27.9"). Data from the EPIC 
cam eras, consisting of one PN back illuminated CCD detec- 
tor (Strudere tai] 120011) and two MOS front illuminated CCD 
detectors ( Tur ner et al.l '200 1 ). have been reprocessed using the 
tasks EMPROC and EPPROC in the Science Analysis Software 
SAS version 9.0. For our purposes it is necessary to rigorously 
select the Good Time Intervals (GTI), representing the periods 
in which the internal background of the cameras was relatively 
quiescent. For this reason we used the SAS task ESPFILT to 
screen data with a high particle-induced background. This task 
was originally de veloped inside the Ex tended Sources Analysis 
Software (ES AS, ISnowden et al.ll2004h and then converted into 
a specific procedure in the SAS. This task fits a Gaussian peak 
to the distribution of count rates, and creates a GTI for those 
time intervals with count rates within the thresholds, defined to 
be +1.5cr from the mean count rate. The GTI filtering has been 
done separately for all the three EPIC cameras. Figure [1] shows 
the results for the screening done on the MOSl data of OBSID 
0202670701, and how the task operates. The upper panel shows 
the histogram for the full Field Of View (FOV) lightcurve in the 
2.5-12 keV band, while in the middle and the lower panels are, 
respectively, the 2.5-8.5 keV lightcurve of the IN FOV region 
and the corner data. The green time intervals are those selected 
by the filtering process and used for the further analysis. The re- 
sults of the screening done on all the selected data sets together 
with the specifics of each observation are reported in Table [1] 

We performed a different GTI selection for OBSIDs 
0402430301, 0402430401 and 0505670101, which were con- 
taminated by very high radiation levels towards the end of the 
respective observations. Indeed, the filtering approach of select- 
ing the count rate inside +1.5cr from the mean count rate can 
lead to a bias in the case of a very strong flare during an ob- 
servation (especially for the PN camera; this does not influ- 
ence the MOS cameras). Moreover when removing the count 
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Fig. 2. Left panel: Fluence map (6.336-6.464 keV) of the Fe-Ko- emission from the filaments in the GC. The image has been smoothed with a 
8-pixel diameter top-hat filter. The grid lines define the Galactic Coordinate frame, where the green circle shows the background region selected 
for the spectral analysis of the clouds marked in the right panel (Section [3.3t . Right panel: A zoom-in on the region of the bright filaments. The 
white (green) ellipses show the regions (sub-regions) selected for the spectral analysis (see TableO; region B has been divided into two subregions 
Bl and B2; region F (GO. 1 1-0. 1 1) into three subregions, Fl (ellipse plus small circle), F2 and F3. The position of Sgr A* is marked by a green star 
at l=-0.056 deg and b=-0.046 deg. 



rates below the mean count rate-1.5cr (see Fig[T]) all the genuine 
low background exposure times in the observation are removed. 
The ESPFILT filtered PN event files for OBSIDs 0402430301, 
0402430401 and 0505670101 have a GTI exposure of 0.3, 4.3 
and 9.2 ks, respectively. We built the 10-12 keV lightcurve of the 
full FOV for these OBSIDs, and selected those intervals with a 
count rate lower than 1.5 counts s ' (for OBSIDs 0402430301 
and 0402430401) and 1.25 counts s"' (OBSID 0505670101). 
The total exposures in the PN cameras of these OBSIDs with 
the new filtering process are 61.5, 48.9 and 74.2 ks, respectively. 

Throughout our analysis we have only selected single and 
double events (PATTERN <4) for PN and up to quadruple events 
(PATTERN<12) for MOSl and MOS2 cameras. For all the in- 
struments a further screening for good events was done selecting 
only events marked as real X-rays (FLAG==0). 

3. Spatial and spectral characterisation of the 
Fe-Ka emission 

3.1. Imaging 

3.1.1. Method 

First we built a fluence (integrated flux over time) map for the 
Fe-Kof line in a very narrow energy window in order to identify 
and locate the brightest Fe fluorescence regions. We proceeded 
by merging all the GTI filtered event files for the three EPIC 
cameras (where available) for every OBSID selected (see Table 
[TJ, using the SAS task EMOSAIC. In building the image of the 
Fe-Kof emission, we assumed a value of E/AE=50 for the spec- 
tral resolution (FWHM) of both PN and MOS1&2 at 6.4 keV, 
corresponding to a bandpass for the Fe fluorescence signal of 
6336-6464 eV. 

An important issue for our analysis is the subtraction of in- 
strument background and, where appropriate, the continuum un- 
derlying the line emission. To quantify the percentage of counts 
in the narrow band 6336-6464 eV due to both components, we 



selected the energy range 3904-6208 eV in the total event file 
and produced 18 adjacent images of the whole FOV, each 128 
eV broad. The spectral window thereby encompassed is free 
of contamination by strong emissions lines, including both the 
Fe-Ko' emission from neutral or near-neutral ions and the Fe-K 
emission from high ionisation lines which characterise the GC 
high-energy thermal components. We then fit the distribution of 
counts as a function of the energy with a powerlaw of the form 
E and further used this best fit to infer the number of counts 
in the continuum for the images in the ranges 6208-6336 eV 
and 6336-6464 eV. Finally, we used the 6272 eV centered image 
properly scaled as background for the Fe-Ko' map. 

3.1.2. Results 

The background- and continuum-subtracted narrow -band (6336- 
6464 eV) image produced by the method outlined in the pre- 
vious section, in effect, represents a Fe-Ko' line fluence map 
(with the provision that the effective exposure varies somewhat 
over the field). This image is shown in Fig|2l The morphology 
of the Fe Ka emission evident in this image is consistent with 
that reported earher by other X -ray satelhtes including A SCA 
dKoyamaet al."1996'), Chandra dYusef-Zadeh et al.ll2007a) . and 
Suzaku (Kovama et al. 200^. The main differences arise be- 
cause of instrumental effects (different angular resolution and 
sensitivity in the hard X-ray domain) or different exposure times. 
Our resul ts are also in good agreement with those reported re- 
cently bv iPonti et al.l ( 1201 Oh . based on an independent analysis 
of a broadly similar set of XMM-Newton observations. 

The maps clearly show an asymmetry in the distribution of 
the Fe-Ko' emission with respect to the position of Sgr A*, with 
all the bright filaments located to the east of Sgr A*, at nega- 
tive Galactic latitudes. On the western side of the GC region, 
the closest molecular complex bright in Fe fluorescent lines is 
Sgr C, located at a projected distance from Sgr A* of approx- 
imately 70 pc. Looking at the radio maps of the GC region 
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Table 2. The central position, sizes (radius or semi-axis) and projected distance from Sgr A* of the regions and sub-regions se- 
lected for the spectral analysis - see Fig|2] Note that region F (GO. 11-0. 11) and also sub-region Fl are comprised of two separate 
components. A cross-reference to the cloud designation employed by Ponti et al. (2010) is also provided. 



Region 


Ponti 


R.A.(J2000) 


DEC.(J2000) 


size (arcmin) 


dpioj (pc) 


The 9 clouds 












A 


MCI 


17:45:52.398 


-28:56:38.15 


0.65x1.08 


10.9 


Bl 


MC2 


17:45:59.592 


-28:57:10.09 


0.91x0.91 


12.6 


B2 


1&2* 


17:46:04.394 


-28:55:56.74 


0.90x0.68 


16.3 


C 


new 


17:45:57.026 


-28:55:22.91 


0.8x0.8 


14.7 


D 


3&4* 


17:46:04.252 


-28:54:39.64 


0.75x0.75 


18.3 


E 


5* 


17:46:10.997 


-28:54:09.50 


1.12x0.83 


21.6 


F 


GO.11-0.11 


17:46:21.627 


-28:53:11.36 


2.75x1.22 


27.2 






17:46:22.615 


-28:55:00.14 


1.68x0.95 




DSl 


new 


17:45:51.294 


-28:53:35.76 


0.77x1.52 


17.0 


DS2 


new 


17:45:43.377 


-28:55:14.99 


1.78x1.16 


12.3 


The 3 F subregions 












Fl 




17:46:17.084 


-28:52:49.45 


1.49x1.01 








17:46:22.443 


-28:52:26.99 


0.59x0.59 




F2 




17:46:28.297 


-28:53:50.10 


1.32x1.01 




F3 




17:46:22.615 


-28:55:00.14 


1.68x0.95 




Large regions 












EDE 




17:45:53.128 


-28:52:26.19 


10.95x7.09 




EDW 




17:45:13.233 


-29:06:46.85 


10.95x7.09 





Notes. The clouds marked with a star are subregions of the Bridge feature identified bv lPonti et alj ( 1201 Oh 



in CO (iQka et alJ|L998!) and CS feboi et aljll999|) rotational 
lines, an asymmetry in the distribution of molecular matter is 
very noticeable, with more concentrated to the east of Sgr A*. 
Nevertheless, a certain amount of cold matter not shining in the 
6.4-keV line is present along the whole plane of the Galaxy, and 
this should be taken into account in explaining the origin of the 
Fe-Ka emission in the filaments between Sgr A* and the GC 
feature known as the Radio Arc. Looking at the left panel in 
Fig|2] the presence of low surface brightness 6.4-keV line emis- 
sion is evident across the whole region between the Radio Arc 
(1»0.2°) and Sgr A* (coded as diffuse blue in Fig|2]i. In con- 
trast, this low-level emission is largely absent at negative longi- 
tudes. We emphasise that this is not an image artifact. In fact, 
the 6.4-keV spatial distribution at faint levels coincides remark- 
ably well with the contours of T eV y-ray emission measured 
by HESS dAharonian et al."2006!), a connection first noted by 
lYusef-Za deh et al. (2007a) using the 6.4-keV line equivalent- 
width image measured by Chandra. Thanks to the much higher 
effective area of XMM-Newton, we are now able to establish that 
this strong correlation apparently extends to regions outside of 
those occupied by very dense molecular clouds (see Section|6|l. 



The bright region in Fig H] located at 1-0.12°, b~0.05° is 
the Arches Cluster; for the study of the 6.4-keV line emission 
fro m the MCs in this region we cross refer to the recent work 
bv lCapelli et al.l (|201 lb) . We have excluded from our study the 
bright feature at l~-0.05°. This is the Sgr A complex, contain- 
ing both Sgr A* itself and the Sgr A East Supernova Remnant 
(SNR). In Table|2]we report the coordinates and the sizes of the 
selected regions of elliptical shape upon which our spectral anal- 
ysis is based. We selected these sky regions solely on the basis of 
Fe-Ka morphology without reference to potential radio molecu- 
lar counterparts (that could help to identify spatially and physi- 
cally connected regions). In doing so, we neglected the possibil- 
ity that some of the filaments which are bright in the Fe-Ka line 
could be part of a single, more extended, molecular complex. 



3.2. Spectroscopy: studying the Fe-Ka variability 
3.2.1. Method 

We carried out a spectroscopic analysis of the Fe-Ka bright 
filaments shown in Fig|2] with the objective of measuring the 
photon flux in the 6.4 keV line and searching for any varia- 
tions in this quantity over an eight-year timeframe. The molec- 
ular complexes bright in the Fe fluorescent line are embedded 
in strong X-ray emission due to the presence along the line 
of sight of unresolved thermal X-ray emission. This thermal 
emission, which permeates the whole GC region can be charac- 
terised as a two-temp erature plasma w ith typically kTia;l keV 
and kT2~6.5 keV dKovama et all 1200 7). For our purposes, ev- 
erything but the 6.4-keV Une emission is considered as back- 
ground radiation. Unfortunately due to the relatively low sur- 
face brightness of the fluorescent emission and the spatial vari- 
ations inherent in the foreground emission, the subtraction of a 
"local background" introduces many systematics in the resul- 
tant spectrum. Therefore, we adopted the approach of modelling 
the Non X-ray Background (NXB) rather than subtracting it. In 
Appendix|A]we describe the different ways of dealing with NXB 
in an XMM-Newton dataset, and specify the methodology we ac- 
tually employed. One consequence of the use of the background 
modelling technique was that for this part of analysis we were 
restricted to the use of the MOS data only. 

We built response and ancillary files with the SAS tasks 
RMFGEN vl.55.1 and ARFGEN vl.76.4. No point source ex- 
cision has been applied. The need for a systematic analysis 
with spectra encompassing very different exposure times and net 
counts, led us to use the Cash statistic (Cash 1979) rather than 
We grouped the channels in each spectrum with the GRPPHA 
tool to ensure at least 1 count/bin. 

We fitted all the spectra with a model which accounts for all 
of the source and instrumental components. The former includes 
the two thermal components noted above, a power-law contin- 
uum plus an associated 6.4-keV Fe line (representing reflection 
and fluorescence in the clouds). We use a second power-law 



4 



R. Capelli et al.: The X-ray lightcurve of Sgr A* over the past 150 years 



Region 


reuUz 


MaiU4 


t>epU4 


Apr07 


MarUs 


AprU9 


A 


2.5+0.3 


2.2+0.2 


2.1+0.1 


2.2+0.2 


2.2+0.2 


2.2+0.2 


Bl 


2.3+0.3 


1.7+0.2 


1.5+0.1 


1.6+0.1 


1.3+0.2 


1.2+0.1 


B2 


0.9+0.2 


0.7+0.2 


1.1+0.1 


2.1+0.2 


2.1+0.2 


2.1+0.2 


C 


0.8+0.3 


0.7+0.2 


0.9+0.2 


1.2+0.2 


1.1+0.3 


1.2+0.2 


D 


0.5+0.2 


0.7+0.2 


0.8+0.1 


1.0+0.1 


1.1+0.2 


2.1+0.2 


E 


2.4+0.4 


2.2+0.3 


1.8+0.2 


2.1+0.2 


2.2+0.3 


2.1+0.2 


r 


1 J. i 3: 1 .u 


i 1 .Z.I!lU.O 




1W.U3:U. J 


Q d+n 7 


Q 1 +n 

y. i 3:U. J 


Fl 


4.4+0.6 


4.0+0.5 


4.4+0.3 


3.4+0.3 


3.0+04 


3.4+0.3 


F2 


4.1+0.6 


3.1+0.4 


3.6+0.3 


3.7+0.3 


3.3+04 


2.6+0.3 


F3 


5.2+0.6 


4.2+0.5 


3.8+0.3 


3.5+0.3 


3.4+04 


3.4+0.3 


DSl 


1.5+0.3 


1.5+0.3 


1.4+0.1 


1.4+0.1 


1.4+0.2 


1.4+0.2 


DS2 


2.8+0.4 


2.3+0.3 


2.4+0.2 


2.1+0.2 


2.4+0.2 


2.5+0.2 



Table 3. Fluxes of the Fe-Ko- line in units of lO^^photons cm"^ s"'. The values are the weighted mean of the M0S1&2 measured fluxes. The 
different columns refer to different datasets (see text). The Apr07, Mar08 and Apr 09 fluxes for the regions D, D S 1 and DS2 are based on MOS 2 
data only, as a consequence of the damaged sustained by CCD6 in MOS 1 on 9 March 2005 jAbbev et al.ll200q) . 



component to model a defined contribution from the instrument 
background, which is not subject to the instrument response (see 
Appendix A). F or a more detailed de scription of these spectral 
components see lCapelli et al. I (l2011bh . Note that in this analysis 
we fixed the temperature of the hot plasma in every spectrum to 
6.5 keV (Kovamaet al. 2007) since this is critical in obtaining 
a consistent measurement of the Fe-Ka line flux across all the 
spectra (different instruments/epochs). 

Throughout our spectral analysis we have considered only 
observations with an exposure longer than 30 ks. If different 
datasets were part of a set of observations carried out over 
a couple of days, the associated spectra were added together 
with the ftool MATHPHA, with the same done for the re- 
sponse files using the ftools ADDARF and ADDRMF. Errors 
have been considered at the 90% confidence level. The result- 
ing temporal sampling was as follows: the OBSIDs 01 1 1350101 
(Feb02), 0202670501 + 0202670601 (Mar04), 0202670701 -t- 
0202670801 (Sep04), 0402430301 + 402430401 + 402430701 
(Apr07), 0505670101 (Mar08), 054750401 + 0554750501 + 
0554750601 (Apr09). A typical result from this spectral mod- 
elling process is illustrated in Fig|3] 

3.2.2. Results 

Our goal is to study the temporal behaviour of the Fe fluores- 
cence in the filaments with a view to constraining the nature 
and location of the energising source. For example, a clear spa- 
tial pattern to the variations might, in principle, help us iden- 
tify a preferred direction for the incoming photons/particles. 
The 6.4-keV line fluxes derived from our spectral analysis for 
the selected regions in the different observing periods are pre- 
sented in Table [3] Here the measurements are calculated as the 
weighted mean of the MOS1&2 values with the errors repre- 
senting the 90% confidence range. Fig|4] shows the resulting 
Fe-Kof flux lightcurves for the different regions covering the 
period 2002-2009. In each of the selected regions the surface 
brightness of the 6.4-keV line lies in the range 5-10x10-*' pho- 
tons cm-- S-' arcmin-^; from the lightcurves it is clear that some 
regions show evidence of substantial variability, while some oth- 
ers appear to remain roughly constant. To quantify this, we fitted 
the 6.4-keV line fluxes as a linear function of time t: 

F64 = A H- B • t, 

where B>0 indicates an increasing Fe-Ka flux (F64). Table |4] 
reports the results for the selected regions and in Fig|4]we show 



data and folded model 




Energy (keV) 



Fig. 3. Upper panel: The spectrum derived for Region B 1 from the 
September 2004 MOSl dataset together with the best-fitting spectral 
model. The lower solid line shows the NXB power-law component. 
The two curves show the net X-ray source spectrum and the source 
plus NXB best fit model (top curve). Three Fe lines are labelled; neu- 
tral Fe-Kor and K/3 lines at 6.4 and 7.05 keV, and the Fe XXV Ko- line 
at 6.7 keV. Lower panel: Ratio between the spectral points and the best 
fit model. 

the best fit linear functions overplotted on the lightcurves, where 
the error bars represent (for this case only) the Icr uncertainty 
as used for the calculation of the ,. From this it is clear that 

^ red 

the temporal behaviour of the Fe-Ko- flux differs markedly from 
filament to filament. In brief, the different regions exhibit the 
following behaviour: 

Region A: the Fe-Ko' flux measured in this molecular complex 
does not show any significant variability on a timescale of 8 
years. The weighted mean flux across the full set of observations 
is 2.16+0.07x10-^ photons cm-^ s-', which translates to a sur- 
face brightness of 9.9+0.03x10-^ photons cm-^ s-' arcmin-^. 

Region Bl: we have measured a decrease of the 6.4-keV flux 
from the sub-region Bl at a confidence of 4.3cr, although a 
linear decrease does not provide a satisfactorily description of 
the light curve (xl^d ~ 1-46, see Table |4|i. The rate of de- 
crease of the Fe-Ko- line flux in this cloud is 8.1+3.1x10-^ pho- 
tons cm-^ S-' year-'. 

Region B2: we have measured a strong increase of the Fe Ka 
flux from the region B2 with a 11.9cr confidence, although the 
high value of the (8.1) for the fit suggests that this rise may 
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Fig. 4. Lightcurves of the Fe-Ka flux from the filaments in the GC. The flux is plotted in units of 10"^ photons cm"^ s"' . The green line represents 
the best fit linear function to the data, whose parameters are reported in Table|4] In each panel we present the significance of the measured gradient 
(B/o-fi) together with thex%j for the best fit linear function. 



not be truly linear. The measured rate of increase of the Fe-Kc 
line flux is 2.4+0.3x10"'' photons cm"^ s"' year"'. 

Region C: the Fe-Ka flux from region C has been measured 
to have an increase at 4.3cr confidence with the overall x^^^ 
for the best fit linear function being 1.2. The rate of increase 
of the Fe-Ka line flux is 7+4x10"^ photons cm"^ s"' year"'. 
Looking at the lightcurve in greater detail, we notice that there 
is a step-wise change at roughly the mid-way point, with the 



weighted means for two segments being, respectively, 0.8+0.1 
and 1.2+0.1x10""'' photons cm"^ s"'. Region C encompasses the 
position of a known X-ray transient, XMMU J 174554.4-285456, 
(iPorquet et al.ll2005h . located at the north-west edge of the circle 
region. This spatial coincidence and the possibility that the tran- 
sient is the illumination source for the Region C cloud will be 
discussed later (§7.1.1). 
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Table 4. Results for the study of the Fe-Ka variability in the selected 
regions and sub-regions (first column). The second and third columns 
give the value of the gradient obtained for the best-fit linear function 
together with its uncertainty (in units of 10"" photons cm"' s"' year"'). 
The fourth and the fifth columns report, respectively, the number of 
sigma for the deviation from the constant behaviour of the best-fit func- 
tion and the reduced for the linear fit. 
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Region D: we have measured a rapid increase of the 6.4-keV 
flux from region D with a confidence of 9.5 cr. The rate of 
growth of the Fe-Kc line flux in this region is 1.6+0.3x10"^ 
photons cm"- s"' year"'. In the lightcurve of region D three of 
the six points are not fitted by the best fit function (see FigH] 
fourth panel); this translates into a high x^red value (9.3, Table 
Ell suggesting the behaviour is not well represented by a simple 
linear increase. 

Region E: the Fe-Kor flux measured in this molecular complex is 
rather constant over a period of 8 years; given the high x^^^ of 
the best fit for the lightcurve of this region {x^,^j-2.9), we cannot 
exclude minor excursions of the flux. The weighted mean flux 
over the set of observations is 2.1+0.1x10"^ photons cm"^ s"', 
which translates to a surface brightness of 7.24+0.03x10"'' pho- 
tons cm"^ s"' arcmin"^. 

Region F: this is the region for which the measurements are most 
precise. The Fe-Ka flux shows a clear decrease with a 7.3 cr con- 
fidence; \htx].ed v^l'is of the best fit function is 0.3, strongly con- 
firming the measured decrease and constraining the functional 
shape of this decrease to be linear The rate of decrease of the 
Fe-Ka line flux is 5.1 + 1.2x10"'' photons cm"^ s"' year"'. To 
improve the spatial resolution, we have divided this region into 
three sub-regions and studied these separately. 

Sub-regions F1-F2-F3: we have measured a decrease of the 
6.4-keV line flux from all three sub-regions with a confi- 
dence of 4.8, 3.4 and 4.2 cr respectively. The coiTesponding 
rates of decrease are 2.1+0.7, 1.5+0.7 and 1.7+0.7x10"'' pho- 
tons cm"^ s"' year"'. The x]-eci values of the best fit linear 
functions are, respectively, 1.3, 3.6 and 1.6; confirming the de- 
crease is a linear function of time. Moreover we point out that 
the rates of decrease from the three sub-regions are compatible 
with each other The agreement between the temporal behaviour 
of both the whole complex and the individual subregions sug- 
gests a common origin for the Fe-Ka emission within a single 
large molecular structure , i.e. the GO. 1 1 -0.11 molecular cloud 
jYusef-Zadeh et alj|2002l) . iTsuboi et al.) (Il997h have show that 
this molecular complex has a shell-like structure that may be in- 
teracting with the western vertical filaments of the Radio Arc. 
Thus, GO. 11 -0.11 could host an extended high energy particle 
accelerator/decelerator. 



Fig. 5. Colour-coded sketch of the molecular clouds and Fe-Ko' 
line flux variability studied in this work. The circles and ellipses 
represent the regions selected for the timing and spectral studies 
(see also the right panel of FiglJ]). The colours are blue (increase 
of the 6.4-keV line flux), red (constant) and yellow (decrease). 
A Galactic coordinate grid, and the positions of Sgr A* (green 
star) and the Arches cluster (cyan stars) are also indicated. 



Regions DSl and DSl: both the low surface brightness regions 
DSl and DS2 show a constant 6.4-keV line signal over a period 
of 8 years. The weighted mean flux within this set of observa- 
tions is 1.42+0.07 and 2.33+0.09x10"^ photons cm"^ s"' re- 
spectively, which translates into a surface brightness of 3.8+0.2 
and 3.6+0.1x10"^ photons cm"- s"' arcmin"^. The x^gj values 
make us confident of the constant nature of the light curves al- 
though, as for region E, we found some little excursions in the 
lightcurve of the region DS2. 

To summarize, we have independently established a pattern 
of Fe-K cy variability similar to what akeady found by Ponti et d] 
(1201 Oh for the clouds A, B2, D, E and F For cloud Bl, our anal- 
ysis implies a decrease of the line flux at a confidence level of 
2.6cr, whereas [Ponti et al] ([2010) report a constant line flux for 
this cloud (designated MC2 in their study); we attribute this sig- 
nificant difference to the use of the background modelling. This 
result has implications for the location of the cloud relative to 
the leading and trailing edges of the ionizing wavefront. We have 
also studied three lower-surface brightness regions (C, DS 1 and 
DS2), the first of which shows an interesting step-wise increase 
in its fluorescent line flux, with the latter two having constant 
light curves. 

In order to investigate the spatial dependency of the observed 
Fe-Ko- variability, and show the complexity of its pattern across 
the inner CMZ, we sketched the distribution of the MCs stud- 
ied in this paper, with a colour-coding which identifies the type 
of variability measured in the 6.4-keV line flux (see Fig|5j- The 
most noticeable increase in the Fe-Ka flux comes from the zone 
bounded by 0.05°<1<0.1° and -0.1°<b<-0.05° (1 and b are, re- 
spectively, galactic longitude and latitude), which encompasses 
the regions B2 and D. Indeed, these molecular filaments show 
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Table 5. Results from the analysis of the PN time-averaged spectra of the Fe-Kc bright filaments. The spectral parameters reported 
are the column density along the line of sight (Nh, in units of 10^^ cm"^), the temperatures (in keV) of the warm and the hot plasma 
components and their normalisations (in units of 10"'^ and 10"'** Jnfn//dV/47rD^ respectively), the column densities intrinsic to the 
MCs (Nhc in units of 10^^ cm"^), the central energy and the intrinsic width cr of the Fe-Ka line (in keV), the 6.4-keV line flux (in 
units of 10"^ photons cm"^ s the EW of the 6.4-keV line (in keV), the excess optical depth t at the Fe-K edge, the power law 
normalisation (in units of 10 photons cm^ s ' keV ' at 1 keV) and the C-stat value and the degree of freedom for the best fit 
model. 
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the clearest evidence for Fe-Kc variability across the selected re- 
gions. Besides these complexes, region C also shows an upwards 
trend in its Fe-Ka surface brightness, wher eas the GO. 1 1-0.11 
cloud shows fading 6.4-keV line emission. iPonti et aP (1201 Ol) 
have suggested that this complex pattern of variability can be 
largely explained in terms of a single flare from Sgr A* illumi- 
nating molecular clouds at differing positions along the line of 
sight. However, as we show below, detailed consideration of the 
spectral properties of the various MCs lead us to conclude that 
this is not the complete story. 



3.3. Spectroscopy: looking for the reflection imprints 
3.3.1. Method 

In order to study the spectral features typical of an XRN, we 
have considered the spectra from the PN camera of the regions 
highlighted in Fig|5](see also Table |2]i. For this analysis we se- 
lected only PN spectra because of the much higher effective 
area in the hard X-ray band (>6 keV). To build the PN spec- 
tra for the Fe-Ka bright filaments, we used the same GTI fil- 
tered event files as employed to construct the fluence map in 
Section 2. We stacked all the spectra (of both the source and the 
background) and the response files using the ftools MATHPHA, 
ADDARF and ADDRMF, neglecting the 6.4-keV line variabil- 
ity. The background region was centred on R. A. = 17:45:42. 289, 
DEC. =-29:07:53. 39 (circle 1 arcmin in radius, see left panel of 
Fig|2|, at the same Galactic latitude as the filaments in order 
to avoid any systematics due to the different intensity of the 
Galactic Ridge emission at different Galactic latitudes. 

The spectral mo del used in the fi tting comprises two ther- 
mal plasmas (APEC, Smith et al. 200 1) with a metallicity fixed 
at 2 times solar (e.g.. lNobukawa et al.l 2010). two gaussian emis- 
sion lines (GAUSS) which account for the Fe-Ka and KJ3 lines 
(the r atio of the line fluxes has been fixed to 0. 1 1, Kovama et all 
l2009l) . and a non-thermal power-law continuum (POW), with a 
slope fixed at r=1.9. The latter is the spectral index obtained 



from a joint fit of the stacked PN spectra of all the MCs using 
the model defined below (where we assume that primary illumi- 
nation source is the same for all the XRN); specifically the best- 
fit value was r=1.9+Q 2. All the previously listed components are 
then subject to a common photoelectr ic absorption by the hne- 
of-sig ht column density, Nh (WABS 1 . iMorrison & McCammonI 
Il983h . 

We model the absorption within the MCs in terms of an in- 
trinsic column density Nh^ (WABS2) plus an additional contri- 
bution to the optical depth of the Fe-K absorption edge at 7.1 
keV. The latter allows the iron abundance relative to the solar 
value (Zpe) in the MC to be greater than 1 (since the WABS pre- 
scription assumes solar metal abundances). In our spectra model 
the extra optical depth t (the free parameter in EDGE, see be- 
low) and Nh(, are tied by the relation: 



T = (Zpe - 1) X O-pe X AFe:H X Nhc 



Nhc 
1022 



X(ZFe- 1)X 0.0116, 



where ApeiH = 3.3 x 10"^ is the (solar abundance) ratio of 
Fe to H atoms and crpe=3.5xl0"20 cm2 is the photoionisation 
cross-section of the Fe-K shell. 

In the event we fixed Zpe at a value of 1 .6 (a choice which 
will be justified in section IJTJ. Note that Nh^ (and r) account 
for the absorption imprinted on the reflected continuum by the 
MC both up to and after the point of scattering. A simple simula- 
tion shows that, to zeroth order, the value of Nh^- so determined 
is comparable to the total column density through the cloud (al- 
though, of course, the exact relation will depend on the geom- 
etry). To summarize, the model used in the spectral analysis of 
the clouds is (in XSPEC format): 

WABSi*[APEC4-APEC+WABS2*(GAUSS+GAUSS-i- 
EDGE*POW)]. 

In the spectral fitting procedure we ignored the channels en- 
compassing the Cu-Ka line (7.8-8.2 keV) in order to minimize 
the effects of the instrumental background on the modelling of 
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The B2 molecular cloud 





Fig. 6. Time-averaged PN spectra of the set of MCs in the present study. Each panel shows the measured data along with the best 
fitting spectral model and (in the lower section) the ratio of these two quantities. The spectral channels encompassing the Cu-Ko- 
line at 8.05 keV (7.8-8.2 keV) have been excluded from the fitting. 



the reflection component. Also, we again employed the Cash 
statistic because of the low signal-to-noise ratio in the spectral 
channels and because its use gave tighter contraints on the de- 
rived parameter values. 

In region F the normalization of the hot thermal component 
was compatible with being zero and therefore we excluded this 
component from the best-fit model. Also, in fitting the spectra 
from regions Bl, C, D and E, we were forced to fix the temper- 
ature of the hot component to kT=6.5 keV (Kovama et al. 2007) 
in order to make the fit converge (see Table [Sj- Although we 
have applied the background subtraction technique in this study, 
all the spectra still exhibit strong residual contributions from the 
two APEC components (as is evident in Fig|6l). 

The specific goal of this analysis was to quantify the intrin- 
sic column density (Nhc) of the MCs and also to establish the 
equivalent width (EW) of the 6.4-keV line with respect to the 
underlying ionising continuum. These important parameters will 
be of use in Sections|4]and|5] 



3.3.2. Results 



ISunvaev & Churazovl(ll998h have shown that in the XRN model, 
the 6.4-keV line emission has to be accompanied by strong ab- 
sorption above 7.1 keV, the minimum energy that a photon must 
have in order to be able to produce a K-shell vacancy in neutral 
Fe via the photoelectric effect. Moreover, because the primary 
source of ionisation is not seen directly by the observer, the EW 
of the line is expected to be very high, i.e. > 1 keV. On the other 
hand, the depth of the Fe absorption edge at 7.1 keV is strongly 
related to the intrinsic column density of the MC, with columns 
in excess of lO^^cm"^ resulting in the imprinting of very strong 
Fe-K edge features on the emergent (i.e., electron-scattered) con- 
tinuum. In Fig|6]and Table|5] we present the results of our search 
for the imprints of reflection on the time averaged spectra of the 
selected GC molecular filaments. We measure a high value for 
the EW of the 6.4-keV line in all the MCs studied. Interestingly 
the largest EWs (~ 1 .5 keV) were those of the B2 and F clouds, 
both of which exhibit pronounced Fe-Ko' flux variability (a de- 
creasing flux in the case of the F cloud and an increasing flux 
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for B2). Another very interesting finding is that the values of 
typically lie in the range 6 to 20x10^^ cm"^, which is sig- 
nificantly h igher than previo us estimates/assumptions for these 
clouds (e.g. lPontiet al.ll2010l) . 

A very important result of our work is the identification of 
the hard non-thermal component associated with the 6.4-keV 
line, which in the XRN model represents the fraction of the in- 
cident continuum Thomson scattered by the molecular material. 
In all the observed spectra we were able to discern an underly- 
ing hard continuum, which via a joint-fit to all the MC spectra, 
was modelled as a power-law with a photon index Fx 1.9. This 
value is very typical of the continuum slopes that characterise 
the spectra of luminous AGN. 

We have also investigated the characteristics of the Fe-Kor 
line, in terms of its energy and width. For this purpose we built 
a time-averaged spectrum for the whole region, stacking all the 
spectra across all the observations. In this stacked spectrum, the 
intrinsic width of the line is 36+3 eV, while its peak is at 6414+2 
eV, slightly higher than the nominal value of 6403 eV for neutral 
iroifl The results from the individual spectra (see Table |5]l show 
the same trend, with the measured Fe-Kc line energy in all cases 
exceeding the theoretical value. Although the EPIC calibration 
is known to be better than 10 eV, the presence of soft protons 
in the FO V might slightly change the energy and the width of 
the lines (iPonti et al.l l2010'). To check whether these are instru- 
mental effects or real features, we measured the same quantities 
for the Cu instrumental line at 8.05 keV (before cutting it out). 
The centroid of the Cu Ka line is at 8047;^2 ^V, while its width 
is 34^^ eV. The line peak is in good agreement with its nominal 
value of 8048 eV. Although a relative shifting between the two 
lines is present (the Cu Ka peak perfectly matches the theoreti- 
cal expectation, whereas the Fe-Ko' is noticeably higher), taking 
into account the systematics and the statistical uncertainties, we 
conclude that no statistically significant line blueshift and broad- 
ening has been measured and that our results are consistent with 
emission from neutral, or close-to-neutral, Fe atoms embedded 
in the material of the cold molecular clouds. The 1 1 eV upwards 
shift of the Fe-Ka line central energy might be the result of very 
modest ionisation of the Fe atoms (Fe 11 to FeX); in fact, the 
ionisation potentials for the low ionisation Fe atoms are between 
few lOeV and few lOOeV, which suggests that these ionisation 
states of Fe might plausibly exist in GC MCs, given the pecu- 
liar environment. We also note that Suzaku measured the Fe-Ko' 
line energy in the molecular filaments to be 6409+1 eV with a 
line width of 33^^ eV, the latter being marginally hig her than the 
expected systematics (^30 eV) dKovama et al.ll20()7h . 

The results of this spectral analysis will be used in Sections 
m and |5] in order to derive a geometrical distribution of these 
nebulae in the GC region and a X-ray lightcurve for the putative 
energising source. 



4. The derived properties of the IVICs 

4.1. Relative Fe abundance in fhe clouds 

In this Section we make use of the measurements of the EW of 
the Fe-Ka line from the spectral analysis of the stacked FN spec- 
tra (see Table|5]l to infer the mean Fe abundance across the MCs 
which compose our sample. The underlying assumption in this 
section is that the Fe fluorescence arises as a result of the X-ray 
illumination of the clouds by one or more X-ray outbursts in Sgr 



' The nominal values for the energy peak of the Fe-Kor and Cu Ka 
lines have been taken from http://physics.nist.gov 




■" 6 



Fig. 7. The cloud-source geometry for the two cloud positions 
{0,n - 0) which produce the same eff'ective scattering. The total 
distance from the cloud to the source (d), the projected distance 
on the plane of the sky (dproj), and the line of sight displacement 
with respect to the Sgr A* plane (dios) are all indicated, as are 
the scattering angles 6 and n - 6. 



A*. To proceed we first require an expression for Fe-line equiv- 
alent width, in terms of the intensity of t he X-ray illumination 
and th e cloud properties. As discussed inl Sunvaev & Churazovl 
( Il998h . the Fe-Ko' line flux (F64) emanating from the illuminated 
cloud can be written as: 



"64 



Q. 
47rD2 



Zpe T"T Ig 



photons cm ^s 



(1) 



In this equation Q is the solid angle (in units of 47r) subtended 
by the diffuse cloud from the perspective of the radiation source, 
D is the distance to the GC, Zpe is the Fe abundance relative to 
solar within the MC, tt is the optical depth due to Thomson- 
scattering through the cloud and Ig is the photon flux from the 
source at 8 keV (in units of photon s"' keV"'). 

In addition the scattered continuum at 6.4 keV measured at 
an angle 6 with respect to the direction of travel of the incident 
radiation is: 



§64 - 



Q. 
47rD2 



Tj (1 +cos^(0))l8 1.172 photons cm^^s^'keV"', (2) 



We can therefore write the EW of the Fe-Ka line with respect 
to the reflected continuum as: 



EW * 



850 Zpe 

[(1 + cos2(6»)] 



eV 



(3) 



In these calculations we assume that the cloud is not optically 
thick along any viewing direction. We also assume the photon 
index of the incident radiation is Fh, ~ 2; the impact of using 
Fin = 1.5 is to increase the EW estimate by ~ 20%. 

If we make the assumption that Zpe is constant across the 
MCs considered in this study, then we can use equation|3]to de- 
termine the value of this parameter. For this purpose we need 
to make a geometrical assumption, namely that the clouds are 
distributed such as to give a reasonably uniform distribution of 
scattering angles within the range 25° < 9 <155°. This is equiv- 
alent to assuming the distance of each cloud along the line of 
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Fig. 8. Left panel: Relative Fe abundance Zpe as a function of the scattering angle 9 for the different MCs obtained by using the 
measured EW and equation (3). The horizontal dashed-dotted line represents the average metallicity (ZFe=1.6) inferred for set of 
MCs. Right panel: Line-of-sight displacement plotted against the projected distance dproj for each MC. The dies value for the DSl 
cloud is a lower limit (dios >9.3 pc). The horizontal dashed line represents the plane of Sgr A* at dio,s=0. The dashed-dotted line 
shows the line dproj=2xdios. We make the assumption that the clouds are located in the region between the dashed and dashed-dotted 
lines, (see text). The permitted dios range (in pc) for each cloud is given in the top-right box. Note: the method does not distinguish 
between whether an individual cloud lies in front or behind the plane of Sgr A* (i.e., at -Hve or -ve dios). 



sight, measured relative to the plane of Sgr A*, is less than twice 
its projected distance from Sgr A*. On this basis, we obtain 
<l-i-cos^(0)>=1.33. The mean Fe-Ko' EW based on the measure- 
ments for the nine clouds reported in Table |5] is 1 .0+0. 1 keV, 
where we have given each measurement a weight proportional 
to the inverse of the error squared. Substituting this weighted 
mean for the EW and the value for <l+cos^(0)> noted above 
into equation[3] we obtain an estimate for the relative abundance 
of Fe in the clouds of ZFe=1.6+0.1 Z©. Note that if we relax the 
geometrical constraint to allow the clouds to be displaced along 
the line of sight relative to the plane of Sgr A* by up to three 
times their projected distance, then the relative abundance esti- 
mate increases by less than 5%. This important result confirms 
a higher than solar metallicity for the MCs and presumably the 
interstellar medium generally in the GC region. 

The geometry of the illumination process is shown in Fig|71 
which also defines the distance along the line of sight behind 
(or in front) the plane of Sgr A* (dios) and the distance of the 
cloud as projected in the sky (dproj), as employed in our work. 
In Fig|7] two clouds are shown at the two angles 9 and 180-0 for 
which the scattering effects are the same. Troughout our work, 
we will always consider a positive dios term for the 180-0 case, 
and negative otherwise. 

Another interesting application of equation |3] deals with the 
geometry of the illumination process. In principle, if the Fe 
abundance is tightly constrained, then the accurate measurement 
of the EW of an individual cloud will allow the scattering angle 9 
to be determined and hence the cloud's position in 3-dimensional 
space to be inferred (albeit with an ambiguity as to whether the 
cloud is positioned in front of or behind the plane of Sgr A*). 
The process is illustrated in the left panel of Fig|8] which shows 
how the inferred relative abundance Zpe varies as a function of 
the scattering angle 9 for the set of clouds in our sample. The in- 
tersection of the horizontal line (corresponding to ZFe=1.6) with 
the plotted curve for a particular cloud determines the scattering 
angle for that cloud (with the positional ambiguity arising from 



the fact that the scattering angle can be 9 or 18O°-0). Although 
there is no such intersection for B2 and F clouds, when the er- 
rors on the EW measurements (which are not shown in Figl8]l are 
taken into we can also derive a range of possible angles for these 
two clouds. The error range on all the scattering determinations 
for all the other clouds can similarly be determined. 

The results of the above analysis are plotted in the right panel 
in Fig[8] This figure shows for each MC the line-of-sight dis- 
placement relative to the plane of Sgr A* versus the projected 
distance from Sgr A* (on the plane of the sky) as determined 
from the scattering angle constraints. In the figure, the MCs are 
positioned behind the plane of Sgr A* but in fact the method can- 
not distinguish whether cloud is in front of or behind this plane 
(see above). Fig[8]also shows that our analysis is self-consistent 
in the sense that the geometrical assumption employed when 
calculating the averaged Fe abundance across the set of clouds 
holds true (i.e., all the MCs lie in the region bounded by the plane 
of Sgr A* and the straight line dios ^ 2xdproj). Unfortunately, 
given the current EW measurement errors, the scattering angle 
determinations results in only weak constraints on the locations 
of a few of the clouds. 

In all of these calculations, we have assumed that the MCs 
are illuminated by the same powerful X-ray source, in a pure 
XRN scenario. Clearly if there is a significant contribution to 
the Fe-Ko' line from an alternative process, such as CR bombard- 
ment, then this would impact on these results. However we show 
in Section|6]that at least for brightest Fe-line emitting clouds, X- 
ray illumination is the most likely excitation source of Fe fluo- 
rescence. Potentially, the method outlined above provides a pow- 
erful tool for the investigation of the illumination geometry of 
XRN, particularly when coupled with studies at radio and other 
wavelengths which may serve to resolve the ambiguity in the 
derived geometry. 
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Table 6. Comparison between the cloud density me asurements. 
The fi rst line reports the Nhc values employed by iPonti et aP 
(l2OT0h . whereas the second lists the results from this paper (see 
Table |5]). All values are in units of 10^^ cm"^. 



A 


Bl 


B2 C 


D 


E 


F 


DSl 


DS2 


4.0 


<2.0 


9.0 


9.0 


9.0 


2.0 






18.4 


10.2 


12.3 5.8 


13.2 


9.6 


9.2 


15.5 


14.5 



4.2. Column densities of tine ciouds 

In the XRN scenario, both the column density of Fe atoms within 
the cloud and the distance of the cloud from the illuminating 
source are the crucial parameters, when estimating the X-ray lu- 
minosity required to produce a given Fe-Ka flux. Previous stud- 
ies have often determined the former using estimates of the hy- 
drogen column density for the cloud, N^^, combined with the 
solar Fe:H ratio and an assumed value for Zpe- In turn, the es- 
timates of have often been inferred from intensity maps of 
CO, CS or other molecular tracers of high density material in 
the inner Galaxy. However, the estimation of the cloud Nh^ from 
molecular line observations is notoriously difficult with different 
measurement techniques often giving different results. For ex- 
ample, the GO. 11 -0.11 cloud (our region F) has been the sub- 
ject of numerous studies over the last decade. Using, respec- 
tively, the intensity of the C S a nd the H^^CO^ emi ssion lines, 
lAmo-Baladron et alJ (l2009h and iHanda et all (12006 ) measured 
two extreme values for the of this MC, namely 2x10^^ and 
10^"* cm"^, respectively. 

A better approach is to measure the column density of Fe 
atoms in the cloud directly from the X-ray spectrum. This is es- 
sentially the approach we have taken in the spectral fitting re- 
ported in Section 13.31 where the Fe column density can be ob- 
tained from the quoted Nhc values as: 

Npe - X ApeiH X (Zpe) = Nh^ X 5.3 X 10"^ cm"^ 

In practice, the fact that we were able to set Zpe = 1.6 and the 
spectral index of the underlying power-law continuum (arising 
from reflection) to r=1.9, resulted in a reasonable robust set of 
Nhc measurements. 

The values we derived for our set of 9 clouds are sum- 
marised in Table |6l where a comparison is made with the ear- 
lier estimat es (based on radio molecular line measurements) by 
iPonti et alJ (1201 Q) for the same MC sample. Clearly the X-ray 
measurements lead to substantially higher estimates of the cloud 
column densities, more similar to the values measured within 
the Sgr B2 complex by Nobukawa et al. (201 1). This in turn has 
a significant impact on the inferred geometry and time delays 
within the cloud system (see below). To conclude, we note that 
the XRN parameters we have found for the MCs within 30 pc 
of Sgr A* are in good agree ment with the ones for the Sgr B2 
MC found bv lNobukawaetaL (,201 1) . i.e. Zpe^ 1.3+0.3 Zq and 
r=2.5+0.6. 



5. The X-ray lightcurve of Sgr A* over the last 200 
years 

5.1. Piacing constraints on tiie iiglitcurve of Sgr A* 

If the number of Fe atoms in a particular MC is known, it is pos- 
sible to use the observed Fe-Ka flux to estimate the cloud's po- 
sition relative to the illumination source, assuming that the cloud 



remains optically thin at the line energy . This method has been 
recently employed bv' Ponti et al.l (1201 0*). who assumed that the 
MCs in their study are all illuminated by a past outburst on Sgr 
A* reaching an X-ray luminosity of -10"''^ erg s"'. However, in 
the previous section 04.2b we noted that the Nh^ values used by 
these authors may be severe underestimates, implying the need 
for a substantial revision of the cloud-source geometry. 

Following the same approach as discussed by 
ISunvaev & Churazovl (Il998h . equation [T] can be used to 
derive an expression for the separation, d, of the MC from the 



lonismg source as: 



d^ = 



10^ Lx TT Zpe 

16;r D2 F64 



(4) 



In this formula tj is the optical depth of the MC to Thomson 
scattering and R is the radius of the cloud. The other parameters 
are the X-ray luminosity Lx of the source (in a nominal 2-10 
keV bandwidth), the relative Fe abundance Zpe within the cloud, 
the distance to the GC D and the measured Fe-Ka line flux F64. 
In applying this formula we make the assumption that the cloud 
is sufficiently symmetric such that we can employ the values of 
Tj and R determined by observation along our particular line of 
sight. 

If we assume that the clouds are illuminated by a p ast out- 
burst on Sgr A* with Lx=1.4xl0^'' erg s"' (e.g. Ponti et al.l 
1201 Oh then we can use equation |4] to calculate the d for each 
of our 9 clouds, since we already know all the other quantities. 
Once the value of d is derived, we can determine the line of sight 
displacement dios of each cloud by using the projected separation 
from Sgr A* dproj on the plane of the sky, (see Table |2]i. Again 
there is an ambiguity as to whether the cloud resides in front of 
or behind the plane containing Sgr A*. Employing this method 
we find that the spread along the line of sight of the MCs is 
almost an order of magnitude larger than the spread in the pro- 
jected distances. Moreover, the inferred line of sight distribution 
of the MCs is incompatible with the 3D distribution we derived 
from the EW relation in Section |4TT] The implication is that we 
must tune (i.e. reduce) the luminosity of the ionizing source to 
match the observed Fe-Kc characteristics of the clouds. 

Ideally, rather than assuming a particular Lx for the outburst, 
we should build a model for the X-ray activity of Sgr A* over 
the last few hundreds years, which accounts for all the Fe-Ka 
variations measured in the MCs in our sample. Unfortunately the 
uncertainty in the position of the clouds along our line of sight 
(relative to the plane containing Sgr A*) hampers this endeavour. 

In the general case, when employing a geometrical model 
such as the one shown in Fig|7] the luminosity of the outburst is 
given as: 



Lx = 2.9 X 10 



39 



F64 (dproj/ sin(0))2 



R2 TT Zpe 



ergs 



(5) 



Also, the delay of the reflected signal relative to the direct 
light from Sgr A* is: 



At = i d, 



pioj 



VTTcot^^ - cot 6* 



(6) 



By allowing 6 to vary over a particular range, one can use the 
above equations to produce the track of Lx versus At matching 
the Fe-Ko- line emission and other properties of the particular 
cloud. Fig|9] shows the result of this analysis, where 6 was con- 
strained to vary within the interval 25° - 155° as in Section l4n 
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However, we can also use the results from Section HTTI which 
are summarised in the right-hand panel of Figl8]to add the further 
constraints in Fig|9] represented by the vertical lines. 

In Fig|9]we have split the clouds into subgroups depending 
on whether they exhibit a decreasing flux (clouds B 1 and F, top 
panel), an increasing flux (clouds B2, C and D, middle panel) 
or a constant level (clouds A, E, DSl and DS2, bottom panel). 
In each of the three panels the different colours refer to different 
clouds with the upper and lower limit ranges on Lx shown by 
the thinner curves of the same colour. Note that in deriving these 
plots we use the Fe-Ka fluxes from the analysis of the stacked 
PN spectra {i.e., the plots are based on the average Fe-Ko' flux of 
the cloud). 

The challenge is to see whether there is a restricted pattern 
of variability in Sgr A*, which plausibly matches all of the con- 
straints implied by Fig|9] Specifically we are interested in those 
Lx-delay combinations within each of sub-groups, which pro- 
vide a consistent description of the Sgr A* variability (i.e., where 
the Lx values for the different clouds within a particular delay 
range are compatible with each other). Of course, given the lim- 
itations we should bear in mind that there can be no guarantee 
that any particular solution is unique. 

The Bl and F clouds: 

In the top panel of Fig|9]we can see that the Lx curves for 
the B 1 and F clouds show a reasonable match around an X-ray 
luminosity of ~4xl0-'^ erg s"' and a time delay around ~100 
years. If these clouds are responding to the same decrease in 
the incident Lx, then we might infer that the scattering angle 
for the cloud F is -95° and that of the cloud Bl neai- to 135°. 
This places the cloud F roughly on the same plane as Sgr A* 
and the cloud Bl 15 pc behind it [ie note that this approach can 
in principle resolve the positional ambiguity alluded to earlier). 
For both clouds the measured excursion of the Fe-Ka flux in 
their lightcurves is about factor of two, implying a commensu- 
rate decline in the Lx of Sgr A* over 7 years. As shown in Fig|9l 
the geometry inferred for these two MCs agrees remarkably well 
with the results obtained earlier (see section 4. 1). 

In conclusion, a drop in the X-ray luminosity of Sgr A* by a 
factor of two, from ~ 4 - 2x lO^*^ erg s"' around 100 years ago 
may explain the decaying Fe-Ka line flux currently observed 
from the B 1 and F clouds. 




1 50 1 00 
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Fig. 9. X-ray luminosity as a function of the time delay for the 
different MCs, as inferred from the measured cloud properties. 
We have divided the MCs into subgroups depending on the type 
of Fe-Ka line variability exhibited. Also plotted are the 90% 
confidence ranges on the Lx calculation. The vertical dotted 
lines represent the constraints on derived in Section l4n Where 
only one vertical line is shown for a cloud, it represents the right- 
hand boundary of the allowed range; in these cases the second 
boundary is situated to the left of the Y axis (i.e., >155 deg), 
and therefore not plotted. 



The B2, C and D clouds: 

In the same way we can try to put together a scenario which 
explains the steady increase in the Fe-Ka flux observed from 
clouds B2, C and D clouds, (middle panel of Fig|9l). Here the 
situation is more complex, since the tracks for the three clouds 
cluster together in a fairly tight fashion along a broad extent of 
the time axis. In order to restrict the possibilities we need to 
bring in extra information as follows. For example, we can per- 
haps safely exclude the time delays < 50 years ago, since the 
requisite activity in Sgr A* might well have been observed on 
this timescale in the early radio observations of the GC (with X- 
ray observations in the 1980s certainly excluding delays of less 
than 30 years). Another good argument which reinforces this as- 
sumption is the level of 6.4-keV line emission measured in the 
50 km s"' and 20 km s"' MCs. These clouds are most probably 
closer to Sgr A* than the ones studied in this paper, and should 
at the present time show a significant Fe-Ka line flux, if illumi- 
nated by an Lx ^lO-'^"'"* erg s"' within the last 30 years. 

Employing the same method as for the B 1 and F clouds al- 
lows us to select a rather narrow range of time delay for the B2, 
C and D clouds; this is the intersection of the three areas in be- 



tween the vertical boundaries in the middle panel of Fig|9] The 
lightcurves of the three regions match well with one another and 
with the results found in Section|4T|for a time delay of between 
70 and 82 years, for an X-ray luminosity in the range 10^^-10^^ 
erg s The corresponding angles, and therefore distances be- 
hind the plane of Sgr A*, for these complexes are about 110°, 
120° and 105° respectively, which translates into a distances of 
6, 8.5 and 5 pc respectively. 

We note that the B2 and D clouds may be sub regions of 
a larger molecular complex called the Bridge by Pon ti et akl 
(2010.) . These authors claimed to have discovered a superluminal 
echo in X-ray reflection, due to the propagation of the ionising 
front within a MC located about 60 pc behind the plane of Sgr 
A*. If we assume that the B2 and D clouds are part of one larger 
complex, then to explain the superluminal Fe-Ka line variabil- 
ity we need to locate the clouds behind Sgr A*, specifically at a 
distance greater than about 8 pc (two times the sum of the lin- 
ear sizes of the regions B2 and D). There is therefore a slight 
disagreement between this requirement and the position inferred 
for clouds B2 and D using the combined results from the spectral 
analysis and the EW-geometry relation. 
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Fig. 10. X-ray light curve of Sgr A* over the past 250 years. The green points a t Lx ~10^ ^ erg s"' show the constraints on the Sgr 
A* activity inferred from the Fe-Ka measurements in the Sgr B2 cloud (see text, Inui et al.f 200 9; Terrier et al. 2010). The red points 
show the current quiescence level of Sgr A*, together with the brightest flare ever measured (lPorauetetalJl2003l P03). The blue 
points in the box show the results of this work. 



To conclude, the observed increase of the Fe-Ka line flux 
in the B2, D and C clouds may be well explained by a short 
enhancement of the X-ray luminosity of Sgr A* in the range 
lO-'^-lO""^ erg s"' about 80 years ago. 

The A, E, DSl and DS2 clouds: 

The bottom panel of Fig|9] shows the Lx curves for the con- 
stant clouds. We can clearly see that the different segments are 
spread over a large area of the plot and it is not straightforward 
to pick out one particular pattern of activity in Sgr A* as the 
most likely solution. In particular, the curve for the DS2 cloud 
does not intersect the curves for the A and E clouds. However, 
the measurements do at least indicate an Lx range rather similar 
to that inferred from the variable clouds. 

From the analysis of the geometrical constraints derived in 
Section 14.11 and plotted in Fig|9] we can try to derive a range 
for the time delay (and therefore position) for the clouds; by in- 
tersecting all the different regions, we propose that the constant 
activity at Lx -lO-*^ erg s"' could have been experienced by Sgr 
A* between 130 and 100 years ago. Assuming a time delay of 



120 years ago, this translates into a distance behind the plane of 
Sgr A* for the clouds of 12, 16, 14 and 15 pc, respectively. 

Finally, in the case of the clouds with constant Fe-Ka flux, 
the displacement between the different curves in Fig|9]might be 
explained in terms of contribution of other excitation processes 
to the total Fe fluorescence (see Section|6]l.] 

5.2. The past lightcurve of Sgr A* 

The analysis of the previous section demonstrates very clearly 
that the Fe-Ka bright clouds within 30pc of Sgr A* are not re- 
sponding to the same outburst on Sgr A* that appears to have 
illuminated Sgr B2. On the contrary, the observed Fe-Ko' line 
fluxes and inferred geometry require an X-ray luminosity of Sgr 
A* more typically in the range lO"*^"-^^ e rg s~^ over t he las t 
100-150 years, in very good agreement with lMuno et al ] (l2007h . 
Within this time frame, upwards and downwards trends lasting 
from a few to ten years also seem to have occurred. 

As a natural step forward, we have attempt to build a long 
term light curve of the X-ray activity of Sgr A*, incorporating 
both published information and our new results - see Fig[TOl The 
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first measurement (the green points in FigfTOt is the X-ray lumi- 
nosity required to explain the strong 6.4-keV line emission in the 
Sgr B2 MC. Almost since its discovery (Kovama et al. 1996), 
there has been a general agreement that the fluorescent of Sgr 
B2 requires a rather powerful X-ray flare on Sgr A*, reaching 
a luminosity of a few xlO^^ erg s"' a few hundred years ago. 
The end of this low luminosity AGN activity has recently been 
dated back to 100_^j years ago (Terrier et al. 2010). Accordingly 
to llnui et al] (|2009|) . tiie 6.4-keV line flux from Sgr B2 was con- 
stant for at least six years (the time of the ASCA first measure- 
ment) before the onset of the decreasing trend in 2000. In FigfTOl 
we highlight the drop in the X-ray activity through green points. 

Between the end of the "high state" activity and the present 
time, Sgr A* appears to have entered a low to intermediate state. 
Fig[TO]shows, in a sketch form, the behaviour inferred from the 
measurements of the clouds within 30 pc of Sgr A*. As dis- 
cussed in the previous section, this variability pattern is not un- 
equivocally determined, although the inferred range of the X-ray 
luminosity is on a sounder footing. The pattern of variability de- 
picted in Fig[TOl does allow the possibility that the regions F and 
B 1 (together with those showing a constant Fe-Ko' line emission) 
might be responding to the declining phase of the flare which en- 
ergised Sgr B2, whereas the other XRN might be illuminated by 
a more recent smaller outburst. 

Further evidence for the long-term downwards trend in the 
Sgr A* light curve is provided by the limits on the Fe-Ko' flux 
from the 50 k m MC. Given th e proximity of this cloud to Sgr 
A* (~ 10 pc) lPontiet"al] ( 1201 Ol) calculate that the mean X-ray 
luminosity of Sgr A* in the past 60 years must have been lower 
than SxlO^*^ erg s"'. For comparison and completeness we have 
also plotted in Fig[TO]the typical level of X-ray emission from the 
SMBH (i.e., ~10^^ erg s~') and the brighte st flare ever detected 
(Lx=3.5xl0^5 erg s ^, iPorquet et al.ll2003h . 

In summary, in the framework of the XRN/Sgr A* scenario, 
we have found that the X-ray light curve of Sgr A* shows a clear 
decreasing trend over the last 150 years. This is not a completely 
smooth variation, with evidence of occasional flaring or, more 
precisely, for periods of relative brightening which last seeming 
longer than ~ 5 years. Given that the light-curve is not mono- 
tonic, a much more comprehensive study of the Fe-Ka proper- 
ties of the MC within the central ~ 100 pc will be needed to 
piece together the full recent accretion history of the SMBH at 
the Galactic Center 



6. The CR contribution to the Fe fluorescence 

6.1. Method 

We have measured the 6.4-keV line flux from two extended re- 
gions to the East and West (EDE and EDW) of Sgr A*, in or- 
der to quantify the difference between the low surface bright- 
ness emission at positive and negative galactic longitudes. For 
this purpose we considered the two elliptical regions shown in 
FigHT] the central coordinates and sizes of which can be found 
in the last two rows of Table|2] In the case of the EDE region, all 
the sky regions coincident with the MCs previously studied have 
been excised, as have the 6.4-keV bright knots in the Arches 
cluster region (see Capelli et al.ll2011bl) . and the Arches cluster 

itself (see iCapelli et al 201 lah . 

We have constructed MOS 1 &2 spectra for both the remain- 
ing EDE area and the EDW region by stacking all the data avail- 
able. Here we apply a background modelling technique in fit- 
ting the four resultant spectra using the same spectral model and 
methodology as described in Section [J!2l 




Fig. 11. A colour-coded and contoured map of the 6.4-keV line 
emission in GC region. The contours levels are scaled in arbi- 
trary units. The position of Sgr A* is marked with the green 
star in the centre of the image, whereas the bright spot close to 
it is the Sgr A EAST SNR complex. The regions selected for 
the spectral analysis lie within the two ellipses marked as EDE 
and EDW. However, the bright MCs studied in this paper and in 
Capelli et al. (2011b), which are all located in the EDE region, 
have been blanked out. 



6.2. Results 

The results of this investigation are summarised in Table |7] 
and in FiglT2] which shows the MOS2 spectra of the EDE and 
EDW regions and the corresponding best-fitting models. There 
are significant residuals near the high ionisation Si, S, Ar and 
Ca lines in the 2^ keV region; these residuals may be due 
to a range of temperatures and/or metallicities within the ex- 
tended plasma which are not well matched by the simple two- 
temperature model which has been applied. The Fe-Ka line flux 
is, however, very well determined in all these spectra. The 6.4- 
keV line fluxes for the EDE and EDW regions are 3.44+0.06 
and 1.91±0. 06x10"'* photons cm"^ s"', respectively. The aver- 
age surface brightness calculated from these values are therefore 
1.71+0.03x10"*' and 7.8+0.3x10"^ photons cm"^ s"i arcmin"^ 
respectively. The substantial higher level of the underlying sur- 
face brightness to the east of Sgr A*, compared to the situation 
to the west, is also apparent in FiglTTI 

The fact that the underlying diffuse Fe-Ka emission at posi- 
tive Galactic longitudes is roughly a factor two more intense than 
that measured at the corresponding negative longitudes is an in- 
teresting result. In particular it can be directly compared with 
the very-high energy y-ray map of the GC region obtained with 
the HESS telescope (lAharonian et al.ll2006l) . The TeV emission 
shows a strong enhancement to the East of S gr A*, with a peak 
close to the centre of our EDE region (see also lYusef-Zadeh et al.l 
I2007bh . 

The TeV emission most likely arises from the hadronic 
interactions of CR with the molecular material in the CMZ, 
which prompts the question as to whether the underlying level 
of the 6.4-keV emission to the east of the GC might simi- 
larly be the result of CR bombardment of the MCs (albeit by 
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Fig. 12. Left panel: MOS2 spectrum of the EDE region, with the residuals in the bottom panel. The data have been collected stacking 
all the MOS2 spectra of the region from all the observations. Right panel: same for the EDW region. 



Table 7. Results from the spectral analysis of the regions EDE 
and EDW. The quoted parameter values are the weighted means 
of those obtained from individual MOS 1 and MOS2 spectra. The 
units of N// are lO^^cm"-, the temperatures of the two plasmas 
are quoted in keV, the normalization of the non-thermal com- 
ponent has units of photons cm"^ s ' keV ' at 1 keV and the 
Fe-KcK line flux has units of photons cm"^ s ' . 





EDE 


EDW 




7.5+0.1 


9.4+0.1 




0.91 +0.01 


0.356+0.002 


kT,„„ 


6.4+0.1 


6.1+0.2 


r 


0.78+0.04 


1.0+0.1 


norm 


5.5+0.1x10^3 


1.61+0.01x10^2 


F64 


3.44+0.06x10-"* 


1.91+0.06x10-"' 


EW 


0.37+0.02 


0.077+0.002 



CR particles of very different energy to those giving rise to 
the TeV signal). A further step is then to conjecture that the 
non-varying emission pedestals seen in some of the 6.4-keV 
bright MCs might also be the result of the same CR bombard- 
ment process. Consider, in this context, the B2 and D regions 
(se e Fig|2l i, which apparently compose the region referred to 
bv iPonti et al. (2010) as the Bridge. The first three data points 
(covering the period 2002-2004) in their 6.4-keV light curves 
(see Table [3] and FiglU are consistent with constant levels of 
1.1 ±0.1 and 0.8+0. 1x10"^ photons cm"^ s"' respectively. These 
measurements translate to surface brightnesses of 5.2+0.5 and 
4.5+0.6x10"^ photons cm"^ s"' arcmin"^, which are compati- 
ble with each other, but higher than the mean 6.4-keV diffuse 
surface brightness calculated for regions EDE and EDW. If one 
then compares the 6.4-keV surface brightness of the B2 and D 
pedestals, the EDE region (after excluding bright MCs) and the 
EDW region, the ratio 6.4/2.2/1.0 is obtained. In a similar way, 
if we look at the contours of the TeV emission map (fig. 12a in 
lYusef-Zadeh et al ] l2007bh . we can see that there is a remark- 
able agreement with this 6.4-keV surface brightness scaling. The 
brightest TeV emission is detected in a region roughly matching 
the physic al extent of the B2 and D c louds (innermost contour of 
Fig. 12a in lYusef-Zadeh et al.l2007bh . whereas the EDE region is 
well delimited by the third contour of the TeV surface brightness 



Table 8. Physical parameters for the MCs and other regions in 
the CMZ region which display Fe-Ko' emission. The table lists 
the cloud column density (in units of 10^^ cm"^), the Fe-Ka line 
flux (in units of 10"^ ph cm"^ s"'), and the sky area subtended 
by the MC or feature (in arcmin^). 



Cloud 


Nhc 


F64 


Area 


A 


18.4!- 


2.5+0.1 


2.2 


Bl 


10.2!?- 


2.3+0.3 


2.6 


B2 


12.3!^-« 


2.1+0.1 


1.9 


C 




1.2+0.2 


2.0 


D 




1.2+0.2 


1.8 


E 


9 6+" 


2.0+0.1 


2.9 


F 




13.0+1.0 


14.6 


DSl 


15 5+3' 


1.5+0.1 


3.7 


DS2 


14.5+2.3 


2 6+"- 


6.5 


EDE 


1.5 


34.0+1.0 


202.9 


EDW 


3.4 


19.0+1.0 


243.9 


N 


4.0 


0.39+0.04 


0.5 


S 


2.0 


0.69+0.05 


0.5 


SN 


2.0 


0.46+0.04 


1.0 


DXc« 


2.0 


0.16+0.05 


0.8 


DXx 


2.0 


0.47+0.09 


0.8 


Sgr B2f 


58.0+3.0 


10.0+1.0 


4.5 


50 km/s* 


9.0 


<1.7 


9.2 


G0.162§ 


2.0 


0.2+0.1 


3.1 


G0.174§ 


2.0 


0.7+0.1 


4.5 


0359.43+1" 


10.0 


0.6+0.1 


7.1 


G359.47t+ 


10.0 


0.9+0.1 


6.9 



Notes , t Ijerrier et al.' t'2010'). t iPonti et al.l fcOlOh . § iFukuoka et id] 
( i2009l) . tt iNakaiima et al.. t200^ . The top section of the table is de- 
voted to the MCs studied in this work, whereas the regions in the sec- 
ond panel are those in the vicinity of the Arches cluster considered by 
Cape lli et al .1 ( 120113) . Parameters for the clouds in the lower section are 
taken from the literature (see the notes below). In the case of the MCs 
with a variable Fe-Kcf line flux, we list the minimum value on record. 
Error ranges for the Nh(. values are not available in all the cases, specif- 
ically when the values are estimates or assumed values. 



distribution, and the EDW region by the next two (lower) con- 
tour levels. 
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Fig. 13. Surface brightness of the Fe-Ka line measured from the 
MCs in the CMZ, plotted against their estimated column density, 
Nhc- The labels indicate the identity of each cloud, while the 
lines represent model predictions for CR electron densities of 
100 and 1000 eV/cm^ (lower and upper lines respectively), for 
assumed slopes of the electron energy spectrum of p=2 (dotted 
lines) or p=3 (dashed lines). CR contribution band as described 
in the text. 



We interpret the above as an indication that CR bombard- 
ment may well play a role in exciting X-ray fluorescence in the 
GC region. Specifically this process could potentially give rise 
to the baseline levels {i.e., emission pedestals) noted in several 
of clouds with variable 6.4-keV emission. In addition, for clouds 
lying outside of the regions where the X-ray illumination is most 
intense, CR bombardment might be the dominant source of fluo- 
rescence excitation (see also Capelli et a l. 201 lb). Plausibly this 
might include some of the more diffuse structures apparent in 
the Fe-Ka maps. 

Unfortunately, it is difficult to disentangle from X-ray stud- 
ies alone, what contribution CR bombardment makes to the total 
6.4-keV line emission. Nevertheless, we can investigate this is- 
sue by assuming (i) that CR electrons with energies in the range 
10 keV - 1 GeV are the cause of the underlying X-ray fluores- 
cence seen in the GC and that (ii) the concentration of these par- 
ticles is relatively uniform across the region of the CMZ studied 
in this paper. Table [8] summarises the physical parameters of all 
the GC clouds and structures which are known to be Fe-Ka emit- 
ters plus some upper limits for several other prominent clouds. 
For the variable clouds, we have considered the minimum 6.4- 
keV line flux recorded to date; for the Sgr B2 cloud we consid- 
ered the Fe-Kg flux value as extrapolated for the year 2010 from 
Fig.4 in lTerrier et all (l2010l) . 

In Figim we have plot the Fe-Ka surface brightness of each 
cloud (in photons cm"^ s"' arcmin"^) versus the estimated col- 
umn density. The contribution of CR electrons to the observed 
Fe-Ka surface brightn ess of various GC clouds has been dis- 
cussed by Yusef-Zade h et al. ( 2007al) . Using the same method- 
ology as .Yusef-Zadeh et al.1 m0 ll\ (in particular their Fig. 13 
and equation 10), we obtain the estimates for the CR contribu- 
tion shown in Fig[T3]as straight lines; here, U and p represent the 
energy density of the LECRe population and the photon index of 



the power law energy distribution of the injected electrons (E ''), 
respectively. FiglT3]has three main regions. The clouds with col- 
umn densities higher than ~5xlO^^ cm"^ and Fe-line intensities 
brighter than ~5xlO"'' photons cm"^ s"' arcmin"^, such as Sgr 
B2 and the bright nebulae studied in the present work, are very 
likely excited by X-ray illumination, i.e., are bona fide XRN. At 
similar brightness levels, but somewhat lower Nhc, lie the MCs 
in the vicinity of the Arches cluste r; as we have shown in our 
previous study (Capelli et al.ll201 fbb . the Fe fluorescence within 
these clouds is most likely energised by particles emanating from 
the Arches cluster itself (CR electrons or protons). The third re- 
gion, lying between the predictions for CR electron energy den- 
sities of 100-1000 eV/cm^, could also plausibly corresponded to 
clouds primarily excited by CR bombardment. 

A remarkable result shown by the plot in Fig[T3]is that, when 
considering the constant (minimum) base level of the variable 
Fe-Ka line emission in regions B2 and D, this falls into the CR 
contribution region delimited by the straight lines in FigfT3l this 
strongly suggest that a CR component might be present in all 
lightcurves studied in this work, and it is of a strategic impor- 
tance to carry out studies on the CR population in the innermost 
CMZ in order to definitely disentangle the contribution due to 
the two sources (CR and Sgr A*), in primis for better constrain- 
ing the past activity of the SMBH. 

Due to the high uncertainties in the parameters plotted in 
this Figure, we do not want to rely too much on it in deriving 
all properties of the MCs we considered, nor putting a net sep- 
aration between the two excitation processes (which could be 
present in tandem in some of the clouds). However, we think 
that some conclusions about the CR contribution to the Fe fluo- 
rescence observed in GC MCs can be drawn already at this stage 
(see FigfT3]i: 

- We can safely conclude that the 6.4-keV line emission mea- 
sured in the A to F complexes is far too high to be only due to 
CR bombardment. Detections of the line variability and the 
absorption feature at the Fe-K edge give further confirmation 
of the XRN nature of these MCs, although a component to 
the Fe fluorescent seems to be present. 

- The EDE and EDW diffuse regions are most likely to be il- 
luminated by LECRs, where the surface brightness gap may 
reflect the averaged dishomogeneity in the CR distribution. 

- The Sgr C MC is very different from Sgr B2 and the XRN in 
the immediate vicinity of Sgr A*; the Fe-Ka surface bright- 
eness is about a factor 20 lower than what measured in Sgr 
B2, and a factor of 10 lower than what measured in the A-F 
clouds. Assuming the same illumination scenario as for the 
Sgr B2 cloud, we calculate that the total distance of the Sgr 
C clumps G359-43 and G359.47 should be as high as 500 to 
700 pc, which seems very unlikely. Moreover, all the X-ray 
observations performed on Sgr C measured the same SB of 
the 6.4-keV line over about 10 years; we therefore propose 
that the description of the Fe fluorescence in the Sgr C cloud 
as due to LECRs bombardment as a likely possibility. 

- The complexes GO. 162 and GO. 174 well lie withi n the CR 
band; as already proposed bv lFukuoka et alJ (1200 9*). the 6.4- 
keV line in the former is likely to be due to CR bombard- 
ment, given the relative vicinity to the radio arc. On the other 
hand, these autors argued that GO. 174 is an XRN, whose 
origin is still not confirmed. For what concerns the plot in 
Fig[T3]we can see that while GO. 162 is well within the CR 
band, G0.174's surface brightness is significantly higher and 
other studies should be performed in order to prove the XRN 
nature of this MC. An Fe-Ka EW of 0.9+0.2 keV is not 
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enough to disprove either hypothesis on the Fe fluorescence 
mechanism. 

Of course, the starting premise of a constant CR concentra- 
tion is far from being correct, but a plot as the one shown in 
FigUHl if nothing else, has the merit of illustrating the order of 
magnitude spread in surface brightness across all the GC clouds 
bright in the 6.4-keV line. 

7. Discussion 

7A. X-ray Reflection Nebulae 

The X RN concept was first sugge sted by Sunyaev et d] 



(Il993h 



Kovama et al.l ( 



1996) 



(120001) . 



and later d evelo ped by 
ISunvaev & Churazovl (1199 8) and iMurakami et all 
In this scenario, a powerful transient X-ray source suppUes the 
primary photons needed to produce the fluorescence observed 
from the GC molecular clouds. The whole phenomenology, 
both the line flux and the topology of the Fe-Ka line emission, 
of the Sgr B2 giant molecular cloud is well explained in this 
scenario. In this context Sgr A* must have been bright roughly 
150 years ago (based on the projected separation of Sgr B2 
and Sgr A*) reachi ng an X-ray luminosity of 2-5x10^^ erg s"' 
dTerrier et al.l l2010b . This is about lO*" times lower than the 
Eddington luminosity for a SMBH with the mass of Sgr A*, and 
so is entirely plausible. 

If Sgr A* is indeed the primary source of the ionising pho- 
tons for Sgr B2 then, by implication, the Fe-Ka emission ob- 
served in the more immediate vicinity of Sgr A* mi ght simi- 
larly serve as a tracer of past outbursts in this source. IPonti et al.l 
(12010) have reported the complex nature of the variability seen 
in the 6.4-keV bright filaments with 30 pc of Sgr A* and have 
suggested a scenario in which these clouds are illuminated by 
the same outburst responsible for the Sgr B2 Fe fluorescence. 

This picture has to date been the favoured explanation of the 
6.4-keV line variability seen in the GC MCs. In broad terms our 
study of the MCs in the immediate vicinity of Sgr A* serves also 
as a confirmation of many of the observational details reported 
by Ponti et al. ( 2010). However, our conclusion is quite different, 
namely that the luminosity required to ionise the MCs studied in 
this paper must be at least one order of magnitude lower than 
the outburst which is currently illuminating Sgr B2. Most im- 
portantly, the derived X-ray lightcurve of Sgr A* over the last 
~150 years appears to exhibit occasional episodes of brighten- 
ing, which last at least 5 years, superimposed on a general de- 
creasing trend. 

The main breakthrough that sets our paper apart from ear- 
lier contributions on this topic is the direct measurement of the 
column density and the Fe abundance of the MCs in our sam- 
ple, which are the two fundamental parameters for the study of 
the XRN behaviour Using this approach we were able to study 
the 3D distribution of the MCs in the innermost CMZ in two in- 
dependent ways, which returned consistent results and therefore 
reinforced our findings. 

First, in Section 14.11 we used our measurements of the EW 
of the 6.4-keV line in order to calculate, within a theoretical 
framework, the average Fe abundance within the selected MCs. 
Our results show that the metallicity of the MCs in our study 
is supersolar with ZFe= 1.6+0.1 Z©, a value which is in very 
good agreem ent with previous estimates (e.g. iNobukawa et alj 
I2010ll2()llh . Secondly, in Section |4j2] we determine the column 
densities of the MCs in our sample directly from the fitting of 
their X-ray spectra. Our X-ray derived Nh^ values are consis- 
tently higher than those estimated from CO-CS measurements. 



Of course, measurements based on molecular lines require the 
assumption of a thermal temperature for the molecular gas which 
is not known a priori, and do not take into account abundance 
variations in molecules such as CS. Another possible way to es- 
timate the of a molecular cloud is the comparison with mm 
continuum maps of dust emission; however, here one meets se- 
vere difficulties in picking out relatively small clouds against the 
general continuum brightness of the GC region. Also the gas-to- 
dust ratio in the GC environment may differ from that found in 
other regions of the Galaxy. 

Both the high value of the EW of the 6.4 keV Fe-Ka line 
and the strong suppression of the non-thermal X-ray continuum 
at the Fe-K edge (at 7.1 keV) supports the view that the bright- 
est clouds in our sample are true XRN. Our determination of 
the spectral slope of the ionising continuum to be F s; 1.9 also 
matches the likely spectral form of the continuum X-ray emis- 
sion emanating from Sgr A* , which illuminates the Sgr B2 cloud 
(e.g. lRevnivtsev et al.ll2004 . 

In the context of the activity of Sgr A* in the recent past, we 
have two important new results: 

- In a pure reflection scenario where Sgr A* supplies the pri- 
mary photons at a constant luminosity, one might expect the 
6.4-keV surface brightness of the densest clouds in the in- 
ner 30 pc region to be between 20-50 times higher that of 
Sgr B2, commensurate with the inverse square fall-off in the 
illuminating intensity. In fact there is no evidence for such 
an effect; on the contrary none of the central region clouds 
have 6.4-keV surface brightness values in excess of that on 
Sgr B2 (see Fig. 9). More quantitatively, we show that the Fe 
fluorescence of the MCs in our sample can be explained in 
terms of the X-ray luminosity Lx of Sgr A* reaching levels 
of between lO-'^ and 10""^ erg s"' over the last 150 years. This 
is more than an order of magnitude lower than that required 
to explain the observed fluorescence of the Sgr B2 cloud. 

- Within the last 150 years the X-ray activity of the primary 
source can be characterised as a declining trend upon which 
occasional episodes of brightening, lasting at least 5-years, 
are superimposed. A scenario in which the observed pattern 
of activity might arise is suggested by the very recent dis- 
covery of a dust clou d on its way towards it s likely final ac- 
cretion onto Sgr A* dGillessen et alj|2012h . Unfortunately, 
there is no detailed knowledge of the physical parameters 
of this cloud, and the X-ray luminosity which might ensue 
following this accretion event is hard to predict with any 
certainty. Nevertheless, on timescales of lO's to lOO's of 
years, accretion events triggered by the capture of passing 
clouds would seem to be a very likely occurrence, partic- 
ularly since the amount of interstellar gas in the accretion 
flow onto the SMBH at th e GC is sufficient, i n principle, to 
produce Lx -10"*" erg s ' dCuadra et alj|2005h . 



7.1 .1 . Association of tine X-ray transient XIVIIVIU 
J 174554.4-285456 witli region C 

The X-ray transient source XMMU J 174554.4-285456 was first 
detected with XMM-Newton during an outburst on 2002 October 
3rd (IPorquet et alJl2005h . Repeated observations of the GC re- 
gion performed with Cha ndra later revealed this source to b e 
bright in June-July 2004 (iMuno et alj|200l IWang et al.ll2006h . 
A closer inspection to the XMM-Newton dataset employed in 
this work revealed the presence of this source in the March 2004 
data. So far, the nature of this peculiar X-ray binary has not yet 
been ascertained. The ratio between the X-ray luminosity at the 
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^ XMMU J174554.4 - 2B5456 

Region C 



Sgr A- 

XMM-NEWTON / EPIC - OBSID 01 1 1350301 



Fig. 14. XMM Newton / EPIC col our image of the X-ray 
transient XMMU J174554.4-285456 dPorquet et al.ll200l . The 
colour coding is as follows: red (0.2-2 keV), green (2-5 keV) 
and blue (5-10 keV). The yellow circle on the top of the image 
shows region C. 



outburst and in the quiescence state is about 10 . For this rea- 
son, we decided to investigate whether the measured fluxes of 
the 6.4-keV line from the region C could be due (in part) to the 
flaring activity of this X-ray binary. 

The observed 6.4-keV line flux lightcurve from the re- 
gion C shows a step behaviour (see Fig|4|. Dividing the six 
epochs into two, we obtain weighted means of 0.8 ±0.1 and 
1.2+0.1x10"^ photons cm"^ s"', respectively. A closer inspec- 
tion to the full band X-ray image from observations 0202670501 
and 0202670601 shows the presence of the X-ray transient 
source XMM U J 174554.4-28545 6 within the bounds of region 
C (see Fig. 15. lPorquet et al.ll2005 ). The best spectral model has 
been found to be an absorbed powerlaw, with a r=1.7±0.2 and a 
very high column density (Nhc) of 14.1;^] "xl0^^cm"^. The flar- 
ing X-ray lumino sity of this source was 1.5xl0^^erg s"' in the 
0.5-10 keV band (iPorquet et al. 2005). 

Assuming the enhancement in the Fe-Kff light curve is due 
to the brightening of XMMU J 174554.4-285456, we calculate 
the requisite distance of the binary to be of the order of 2 pc. 
Given the high Nh^ value measured towards this source, we 
favour the geometry in which the MC lies in the line of sight 
between us and XMMU J174554.4-285456. In this setting, the 
measured time delay between the first reported flare on the bi- 
nary (October 2002) and the onset of enhanced 6.4-keV flux (4.5 
years), the size of the cloud (3.7 pc) and the distance of the X-ray 
source (2 pc) represent a consistent picture. 

In summary the Fe fluorescence of region C may, in princi- 
ple, be explained without recourse to the XRN/Sgr A* paradigm. 
Its base (or pedestal) level of Fe-Ka emission might arise from 
CR particle bombardment (see Section |6]l upon which is su- 
perimposed a component due to X-ray irradiation of the cloud 
by a (very) nearby transient source XMMU J174554.4-285456. 
Potentially this is the second XRN posit ively associated with an 
X-ray transient source in the GC region (ICapeUi etalJl20im . 

7.2. CR Particle bombardment 

This m odel has been proposed and discussed by Valinia et alj 
(I2000I). [Predghl et aL (2003), Yusef-Zadeh et al. (2007a) and 
iDogiel et al.l ( l2009l) . In this context, the primary source of the Fe 
fluorescence are CR particles, potentially with a variety of com- 



positions, origins and energies. While the photoionisation cross- 
section is a steep function of the energy of the primary photon, 
the cross-section for collisional ionisation is a relatively smooth 
function of the energy of the incident particle. For electrons with 
energies in the range 10-100 keV and protons with energy be- 
tween 10 and 100 MeV, the cross section for collisional ionisa- 
tion of Fe remains relative close to its maximum value, i.e. these 
are the energy ranges over which the particle bombardment most 
likely occurs. However, the number of particles will fall sharply 
with energy so one might not expect the major contribution to be 
due to the highest energy component. 

Although electrons have been historically favoured as the 
source of ionisation of MCs in the CMZ, electrons with a kinetic 
energy of some tens of keV can barely penetrate the target clouds 
since they are stopped by a column density of only ~ 10^'cm"^. 
As a consequence, fluorescent emission via electron bombard- 
ment most likely takes place within the surface layers of a dense 
cloud. On the other hand, protons have a much larger penetrating 
power than electrons; a 30 MeV proton is stopped by a column 
density higher than lO^^cm"^, a value approaching the column 
densities found in the GC region. Accordingly, protons may be 
as, or more, important than electrons as a source of ionisation 
of MCs, although peculiar locations (the Radio Arc, the Arches 
cluster) might reve al different me c hanism s at work. 

Very recently, ICapelli et al ] (12011b') showed that in the 
Arches cluster region there are MCs similar to those studied in 
this paper in terms of their Fe-Ka surface brightness, size and 
inferred density. In fact, it is very likely that the bulk of the un- 
derlying fluorescence in these clouds is induced by particle bom- 
bardment, with an estimated CR energy density of about 50-100 
eV/cm^, that is two orders of magnitude higher than that esti - 
mated for the Galactic plane in general (see lCapelli et al.l201 Ibl) : 
the Arches cluster is, of course, the most likely site of the requi- 
site particle acceleration. 

Our results confirm that a major contribution to the Fe-Ko- 
line emission in the MCs in our sample comes from X-ray pho- 
toionisation. However, we argue that a significant contribution 
from CR bombardment might be present in the spectra of these 
XRN - as discussed in Section|6] 

To support this hypothesis, we presented in FigfT3] possible 
evidence for the underlying particle ionisation component. In a 
plot of the Fe-Kc surface brightnesses against the cloud column 
density of the MCs, we were able to distinguish diff'erent regions 
of parameter space in which X-ray illumination and CR particle 
bombardment might operate as the primary source of the Fe flu- 
orescence. However, much more information is required in re- 
lation to the CR energy density and spatial distribution in the 
innermost CMZ is required if we are to understand the relative 
importance of the various excitation processses. 

The Fe-Kc line emission from particle bombardment is gen- 
erally expected to show a constant lightcurve. However, it is at 
least conceivable that time variable magnetic fields in the fila- 
ments and threads and/or time variable CR fluxes might induce 
measureable changes over timescales upwards of a few years 
in some compact molecular structures. Evidence for such be- 
haviour could come from variability studies of the radio contin- 
uum in narrow radio emitting filaments confined by magnetic 
pressure, although of course in such cases the particles will be 
electrons with relativistic energies emitting synchrotron radia- 
tioii; 

lOka et akl (1201 0*) discovered high velocity compact clouds 
in the GC region which show a CO J=3-2/J=l-0 intensity ratio 
higher than that typically found in the Galactic disc. These fea- 
tures are likely to be knots of shocked molecular gas, heated by 
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energetic events such as supernovae explosions. The distribution 
of these clumpy structures correlates well with the location of the 
molecular filaments we have studi ed in this work. The same is 
true for the Sgr C molecular cloud (lOka et al.ll2010l) . This corre- 
lation suggests a close relation between these shocked molecular 
complexes and the Fe-Kc line emission, since shoc ks internal to 
the clo uds act as particles accelerators. Moreover, iMuno et aP 
(l2008h showed that there are many physical processes active in 
the GC region which appear to produce diffuse X-ray features. 
For example, these authors discovered 20 candidate pulsar 
wind nebulae within 20 pc of Sgr A*, which are among the best 
established sites of CR acceleration. The particle bombardment 
hypothesis places few direct constraints on the cloud/particle 
flow geometry, given that particle will travel along the local field 
lines. Nevertheless this scenario does provide at least a plausi- 
ble explanation of why the brightest Fe fluorescence regions are 
sometimes located on the farside of the c loud as viewe d from 
Sgr A* (for example, in Sgr C (.Yusef-Zadeh et al.ll2007"al) . 

8. Conclusions 

We have studied the X-ray properties of selected molecular 
clouds and filaments to the east of Sgr A* at a projected distance 
of ^8-30 pc. 

- Fe-Ka topology: we have demonstrated that the Fe-Ka emis- 
sion delineates both compact and more diff'use structures. We 
have studied the temporal and spectral properties of 9 clouds; 
three of which are newly studied (regions C, DSl and DS2). 

- Fe-Ka variability: significant variability is seen from some 
of the molecular clouds, although the pattern is by no means 
a simple one. Our Fe-Ka flux and variability measurements 
agree reasonably well with previous published results, al- 
though in the case of one cloud (Bl) we measured a decrease 
of the Fe-Ka line flux, in contrast with what found previ- 
ously. 

- Fe-Ka surface brightness: we measure the surface brightness 
of the 6.4-keV line to be of the same order of magnitude in 
all the molecular filaments we have studied. Typically the 
observed variability only involves a factor of two change. 
In the Sgr A*/XRN scenario, the surface brightness of the 
filaments might be expected to decline with distance from 
Sgr A*. This issue is resolved by our finding that the X-ray 
luminosity of Sgr A* has declined substantial over the last 
150 years. 

- Fe-Ka EW and Fe-K edge measurements: we measured a 
high value of the EW of the 6.4-keV line from all the MCs in 
our sample, consistent with an origin of the bulk of the flu- 
orescence in photoionisation. Moreover, all the clouds show 
an absorption feature at the Fe-K edge energy of 7. 1 keV. 

- Spectral hardness of the reflected continuum: the spectral 
shape of the reflected/ionising X-ray continuum emission as- 
sociated with the Fe fluorescence is found to be F ^1.9, a 
result compatible with the XRN/Sgr A* scenario, since this 
is a value rather typical of that found in AGN. 

- Fe abundance: we use the measured EW values of the Fe- 
Ka line in our sample of MCs to calculate the average Fe 
abundance (relative to solar). We show that ZFe=1.6, a result 
which is consistent with the general finding that a higher than 
solar metallicity characterises all the ISM phases in the GC 
region. 

- Column density through the MCs: we have employed a rela- 
tively direct method to calculate the hydrogen column den- 
sity, Nhc, within the MCs. This is based on the joint mod- 



elling of the low-energy absorption and the absorption at the 
Fe-K edge imprinted on the incident continuum in its pas- 
sage through the cloud (both up to and after the point of scat- 
tering). This approach benefits from prior knowledge of the 
relatively iron abundance, Zpe, in the cloud. The column den- 
sities so calculated are close to 10^^ cm"^ and, in most cases, 
significantly higher than those inferred in previous studies, 
which made use of radio molecular line measurements. 

- The past X-ray activity of Sgr A*: based on our study of the 
6.4-keV line emission properties of the MCs in the inner 
30 pc of the GC region, we outline a model for the X-ray 
emission of Sgr A* encompassing the last ~150 years. Over 
this period the X-ray luminosity has declined from an appar- 
ent peak of ~10^^ erg s"^ roughly 150 years ago, to 10^^^"^^ 
erg s"' perhaps 100 years ago, down to typically lO^"'"^^ 
erg s ' at the present time. This decreasing long-term trend 
has, however, been punctuated by counter-trend episodes of 
brightening by factors of a few over timescales in excess of 
~ 5 years. 

- Cloud C - an XRN energised by a transient source?: we 
have found that the 6.4-keV line flux variability measured 
from cloud C could also be associated with the transient 
X-ray source XMMU J 174554.4-285456. In this scenario, 
the Fe fluorescence in this complex is composed of a steady 
non-zero level, possibly produced by the interaction of CR 
with the ambient gas plus a superimposed contribution due 
to the localised variable X-ray irradiation. If this picture 
is confirmed, this is the second XRN found in the GC re- 
gion which has been illuminated by a transient X-ray source 
dCapeUi etaDIMTbl) . 

- Low surface brightness 6.4-keV line emission: we have found 
a good correlation between the TeV emission and the Fe- 
Ka line emission observed on arcminute scales in the inner 
CMZ. Specifically we measured diffuse, low surface bright- 
ness emission to the east of Sgr A*, which seems to perme- 
ate the whole region between the Sgr A* and the giant Radio 
Arc. A pedestal {i.e. constant) fluorescence component might 
in fact be present in all the lightcurves which show a vary- 
ing Fe-Ka line flux. We have also shown that this diff'use 
low surface brightness 6.4-keV line emission correlates quite 
well with the contours of the brightest features seen in TeV 
y-rays. This suggests that CRs may contribute to the ionisa- 
tion and fluorescence of cold gas in this region. 

- CR contribution: we presented arguments in favour of a CR 
contribution to the total Fe fluorescence in the GC MCs. In 
particular, LECR may be responsible for the constant of 6.4- 
keV line emission observed in several structures (including 
the Bridge), as well as in other regions of the CMZ, like the 
Arches cluster, the Radio Arc and Sgr C. We suggest that 
further study of the Fe-Ka surface brightness versus Nh(, re- 
lation for GC clouds will help reveal the contribution of CR 
to the Fe-K fluorescence observed in the GC region. 

We are far from unveiling the full mystery of the Fe fluorescence 
observed in the GC region. It will be of great importance in the 
future to continue to monitor the Fe-Ka line emission from the 
MCs in the CMZ so as to better characterise the past activity of 
Sgr A*. Future X-ray observations can also complement stud- 
ies in different energy domains, from the radio band (continuum 
and lines) up to the TeV regime, which aim to fully chart the 
influence of CR particles on the GC environment. 
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Appendix A: Dealing with tlie EPIC detector 
background 

Within an XMM-Newton EPIC dataset, two main background 
components can be defined, namely particle-induced events 
and photon-induced events. Solar-flare particles and CRs are 
responsible for the continuum element of the non X-ray 
Background (NXB), as wel l as for the strong instrumental lines 
jLeccardi & Molendill2008l) . A further minor component within 
the particle-induced background is due to Quiescent Soft Protons 
(QSP). These are particles which enter the optics and are focused 
onto the detector; they are mostly filtered during the GTI selec- 
tion, but is likely that there is a low-level residual component 
in the screened event file. In the case of the photon background 
from the sky, the main components are the hard X-ray emission 
along the plane of the Galaxy (the Galactic Ridge emission both 
in the foreground and extending beyond the GC), and also the 
extragalactic Cosmic X-ray Background (CXB). We did not take 
into account t he Galactic Halo since its e mission is negligible 
above 2.0 ke V (iLeccardi & Molendil200"8h . The difficulty of car- 
rying out a systematic analysis with a careful characterization of 
the background is that the particle-induced background, which is 
the dominant component in the hard X-ray band, changes both 
spatially on the detector and in time throughout the duration of 
the mission. Given the focus of our investigation, there was a 
clear requirement that all the datasets should be analysed in a 
very systematic fashion. This led us to consider the advantages 
and disadvantages of the standard methods of dealing with the 
background in XMM-Newton observations. 

Local Background subtraction: this procedure consists of select- 
ing a region in the same image/dataset close to the source of 
interest; it is certainly the best choice for the analysis of a bright 
point source. The fact that we are dealing with extended sources, 
which are close to each other, forces us to select a region well 
separated from the position of the filaments under study. This 
may, potentially, produce a large discrepancy between the inten- 
sity of the instrumental lines in the source and in the background 
spectra, especially in the PN camera {e.g. at the Cu line at 8.05 
keV), generating strong systematic errors which are difficult to 
deal with. As a further difficulty, the presence of the Galactic 
Ridge Emission throughout the whole FOV prevents us from se- 
lecting a region for the background with a low surface bright- 
ness. 

Blank Sky Fields: Blank sky files are constructed using a super- 
position of pointed observations o f regions of the sky without 
contamination by bright sources dCarter & Re ad 2007). They 
contain most of the background components listed above, like 
the particle induced NXB and the CXB. To study an extended 
source, we tried to cast the blank event file onto the sky at the 
position of our observations using the skycast script. In doing 
so, we are able to select the same detector region for the produc- 
tion of both the source and the background spectra. This method 
has the great advantage of completely avoiding the possibility 
of having different instrumental lines intensities, but the vari- 
ous components of the background can all vary in different ways 
over long time periods. So far, the blank sky files available are 
for observations done before revolution 1200, that is early 2007 
(private communication, A. Read). We tested the public blank 
fields in our spectral analysis. While for revolutions earlier than 
1200 the instrumental background subtraction was feasible, the 
spectra obtained with later data still show strong contamination, 
which is easily noticeable in the intense Cu Ka line at 8.05 keV 
in the PN camera. Since we are using XMM-Newton observa- 
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tions taken over a period of ~10 years, the blank sky event files 
technique cannot be used for our purposes. 

Background modeling: this technique requires a careful char- 
acterization of all the background components. The strong ad- 
vantage of modelling the background rather than subtracting 
it, is that many of the systematics may thereby be avoided. 
IPe Luca & Molendil (|2004) and iLeccardi & Molendil (l2008h 
studied in detail the properties of the particle-induced back- 
ground for the MOS cameras, giving a spectral shape to both 
the continuum and the line emission. Modeling the background, 
instrumental and cosmic components we can, in principle, ac- 
count for all the photons collected by the instrumentation and 
make every photon count. 

On the basis of the above, we decided to follow the back- 
ground modelling technique for measuring the temporal vari- 
ability of the Fe-Ko' emission. With this choice we did not con- 
sider further the EPIC-PN data, but worked only with MOS cam- 
eras. Up to the present, a self-consistent and exhaustive charac- 
terization of the background for EPIC-PN camera has not been 
developed, mainly because the fraction of out-of-time events 
in the PN CCDs is not negligible and the Out FOV region 
is much smaller compared to the one in the MOS detectors. 
Moreover, MOS cameras are better suited to the study of ex- 
tended sources because of the lower level of the instrumental 
background (which is most likely a consequence of the different 
CCD structure). However, for the analysis of the stacked spectra 
the excellent hard response of the PN camera outweighed these 
considerations and for this study we employed the PN spectra 
and local background subtraction. 
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Table 1. Specifications for the selected OBSIDs: MODE/FILTER combination used for each of the pointings, and GTl compared to 
the total exposure for each instrument. F=FuU Frame MODE; E=Extended Full Frame MODE; T=Thick filter; M=Medium Filter. 
OBSID 0506291201 has PN in Timing MODE. 







Instrument specifics 


GTl & Exposure (ks) 


OBsm 


Obs Date 


PN 


MOSl 


MOS2 


PN 


MOSl 


MOS2 




yyyy-mm-dd 


mode/filter 


mode/filter 


mode/filter 


GTl/exp 


GTI/exp 


GTI/exp 


0111350101 


2002-02-26 


F/T 


F/M 


F/M 


38.590/40.030 


42.262/52.105 


41.700/52.120 


0202670501 


2004-03-28 


E/M 


F/M 


F/M 


13.320/101.170 


33.070/107.784 


30.049/108.572 


0202670601 


2004-03-30 


E/M 


F/M 


F/M 


25.680/112.204 


32.841/120.863 


35.390/122.521 


0202670701 


2004-08-31 


F/M 


F/M 


F/M 


59.400/127.470 


80.640/132.469 


84.180/132.502 


0202670801 


2004-09-02 


F/M 


F/M 


F/M 


69.360/130.951 


94.774/132.997 


98.757/133.036 


0402430301 


2007-04-01 


F/M 


F/M 


F/M 


61.465/101.319 


61.002/93.947 


62.987/94.022 


0402430401 


2007-04-03 


F/M 


F/M 


F/M 


48.862/93.594 


40.372/97.566 


41.317/96.461 


0402430701 


2007-03-30 


F/M 


F/M 


F/M 


32.337/32.338 


26.720/33.912 


27.685/33.917 


0505670101 


2008-03-23 


F/M 


F/M 


F/M 


74.216/96.601 


73.662/97.787 


74.027/97.787 


0554750401 


2009-04-01 


F/M 


F/M 


F/M 


30.114/38.034 


32.567/39.614 


33.802/39.619 


0554750501 


2009-04-03 


F/M 


F/M 


F/M 


36.374/42.434 


41.376/44.016 


41.318/44.018 


0554750601 


2009-04-05 


F/M 


F/M 


F/M 


28.697/32.837 


37.076/38.816 


36.840/38.818 



